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§ 1. INTRODUCTION 
1.1. Cooperative Assemblies 


Tue physical phenomena to which statistical mechanics has been applied 
can, broadly speaking, be divided into two categories. In the first category 
the microscopic systems of which the statistical assembly is composed 
can be regarded to a good approximation, as non-interacting. The funda- 
mental results of statistical mechanics relate thermodynamic functions 
to the energy levels of the total assembly, and the latter can be determined 
quite simply from the energy levels of the individual microscopic systems. 
If the systems obey Bose-Einstein or Fermi—Dirac statistics an elementary 
combinatorial problem must be solved, but this involves no serious 
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difficulty. Of course there must be sufficient interaction between the 
systems for equilibrium to be reached during the time considered, and for 
statistical mechanics to be applicable; but we can reasonably assume 
that equilibrium can be achieved by interactions which make a negligible 
direct contribution, to the energy levels of the assembly. 

Examples of phenomena of this kind are the specific heats of gases, 
and solids, chemical reactions in gases and their equilibrium constants, the 
energy distribution, of black-body radiation, paramagnetism, the elemen- 
tary electron, theory of metals, and the condensation, of a perfect Bose— 
Einstein gas. In the case of solids, although the interactions between 
individual atoms or molecules play an essential part, their departures from 
equilibrium are quite small over a substantial temperature range. Hence, 
by a transformation to normal coordinates the solid can be treated as an 
assembly of non-interacting harmonic oscillators. With the exception 
of the Bose-Einstein condensation, the thermodynamic properties of 
assemblies of the first category are smooth and continuous, and the results 
obtained by the application of statistical mechanics to these assemblies 
can, be considered as most satisfactory. 

In regard to assemblies of the second category, however, the situation 
is very different. For the most part the physical phenomena which they 
attempt to account for contain discontinuities corresponding to phase 
transitions of various kinds. Examples are the condensation of gases, the 
melting of solids, the properties of the liquid state (including the phenomena 
associated with the critical point), the properties of solutions and the 
conditions for phase separation, the form of solubility curves, various 
types of A-point transition (order disorder transitions, the liquid helium 
transition, the superconducting transition), and ferromagnetism. For 
assemblies of this kind the interactions between individual microscopic 
systems play an essential part, and cannot be removed by elementary 
transformations. 

The energy levels of the assembly are not simply related to those of the 
microscopic systems, and the interactions cause the cooperation of large 
numbers of systems giving rise to phenomena such as those listed above. 
The terms ‘cooperative phenomena’ and ‘cooperative assemblies’ are 
therefore often used to describe the second category. The exceptional 
case of a Bose-Einstein condensation referred to above in the first category 
can, also be considered as arising from an attractive interaction between 
the particles, Bose-Einstein statistics being equivalent to such an 
attraction. 

1.2. Simplified Models and Approximations 

The mathematical problems associated with cooperative assemblies are 
extremely formidable. To facilitate calculations models have been 
introduced in which the molecular interaction is simplified, but which, it is 
hoped, maintain the cooperative characteristics of the original assembly ; 
examples are the use of hard sphere molecules in elucidating the properties 
of liquids, and the neglect of all but nearest neighbour interactions in 
problems on crystal lattices. 


M2 
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For the most part even the simplified problems still show little hope of 
exact solution, and approximations have been introduced (e.g. the super- 
position approximation in the theory of liquids Kirkwood 1935, Kirkwood 
et al. 1950, Born and Green 1946, 1949). To make effective use of such 
approximations it is very desirable to have some means of estimating the 
error involved, and hence of delineating the region in which the approxi- 
mation is reliable. It has become fairly clear that in the critical region 
of phase transitions single approximations are unreliable, and the only 
hope of deriving information is by assessing the trend of a series of successive 
approximations; but for the theory of liquids even the simple approxi- 
mations involve intricate and complex calculations, and there have been, 
few attempts at higher order approximations. Recently numerical 
computations on an assembly of hard spheres (Alder and Wainwright 
1957, Wood and Jacobson 1957) have given very strong indications of a 
first-order phase transition at a definite specific volume; no such firm 
conclusion, could be drawn from approximate treatments. 

In regard to cooperative phenomena in crystals the situation is appre- 
ciably better. A particular set of exact solutions exists for two-dimensional 
models. These exact solutions are most valuable, and demonstrate clearly 
that even the simplest type of interaction causes the nature of the solution 
to change completely, and gives rise to characteristic singularities in the 
thermodynamic functions of the assembly. For three-dimensional 
models, although exact solutions still seem far off, approximate methods 
and series expansions have been developed to quite a high order in recent 
years. It now seems possible to specify with confidence many of the 
important thermodynamic properties of simple models, and hence, if 
there is disagreement with experimental results, this must be due to the 
nadequacy of the model and not to a particular approximation. 


1.3. Scope of the Present Article 


The present article is concerned with cooperative phenomena in crystals. 
In § 2 the general theory of regular assemblies is developed, and the basic 
mathematical problem which they pose is formulated. The relationship 
or various simplified physical models to this problem is established. 
§ 3 attempts to give an account of the exact information which is available. 
The main emphasis is on the physical conclusions, but some effort has 
also been made to summarize the important mathematical methods. 
§4 deals with approximate methods, particular consideration being paid 
to the systematic development of successive approximations. By care- 
fully examining the trends of these approximations an assessment is made 
of the physical properties of three-dimensional models. The approxima- 
tions considered in § 4 include the use of finite numbers of terms of exact 
series expansions, and technical problems connected with the evaluation 
of terms of such expansions are dealt with in §5. These terms can be 
related to the topological properties of finite configurations on crystal 
lattices, and special methods have been devised in recent years of deter- 
mining these properties. The results of §5 are also of relevance to a 
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detailed discussion of ferromagnetism and antiferromagnetism on the 
Ising and Heisenberg models which the writer hopes to include in a subse- 
quent publication. The choice of material for §4 has of necessity been 
somewhat arbitrary. The general criterion has been to include properties 
of general significance to the theory of phase transitions, and to defer 
properties which are more specifically magnetic (e.g. the magnetic suscep- 
tibility at temperatures above the Curie point) to the subsequent publica- 
tion. The latter will also contain some exact results of relevance to 
magnetism which have been obtained recently, and will endeavour to 
provide a summary of other cooperative problems on crystal lattices 
which are not dealt with in the current text. 

The literature on cooperative phenomena in crystals is very extensive 
and no attempt has been made to give adequate reference to all research 
work on the subject. The main reason for writing the article was the 
feeling that a stage had now been reached at which definite conclusions 
could be drawn regarding the behaviour of simple models of such coopera- 
tive phenomena. The article attempts to explain the major developments 
leading to such conclusions. 

The question of suitable notation has been a continual source of difficulty 
because of the many different symbols which had to be used. A list of the 
principal symbols in order of their introduction is therefore provided. As 
far as possible the individual chapters are self-contained. 

Throughout the preparation of the article the writer has drawn freely 
on the help of his colleagues Dr. M. E. Fisher and Dr. M. F. Sykes, and his 
former research student Dr. J. M. Burley; all of these placed at his disposal 
the results of their investigations before publication. More specifically 
he is indebted to Dr. M. E. Fisher for a careful and critical examination 
of §3, to Dr. D. M. Burley for advice in connection with § 4, and to Dr. M. F. 
Sykes for collaborating in the writing of §5. Work on the article began 
during the tenure of a visiting professorship at the University of Maryland ; 
the facilities granted by this University are also gratefully acknowledged. 
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§ 2. ReautaR ASSEMBLIES 
2.1. Introduction 


For many of the cooperative assemblies mentioned in §1 a classical 
approximation is adequate, and we can consider our molecular constituents 
as point particles with an interaction energy between any pair which 
depends on their distance apart. The total energy of the assembly is 
then the sum of all such possible binary contributions. Clearly the 
properties of the assembly are determined by the form of the intermolecular 
potential energy, and we should have to take account of all possible 
distances of interaction. However, the problem is simplified if the 
molecules are situated on a crystal lattice, since only values of the 
interaction at a discrete series of points need be considered. If further 
we make the assumption that the intermolecular forces are of sufficiently 
short range for all but nearest neighbour interactions to be neglected, we 
need consider the value of the interaction energy at one point only. This 
provides us with the simplest possible model of a cooperative assembly 
to which the term ‘regular assembly’ has been assigned by Fowler and 
Guggenheim (1939). 

A regular assembly thus possesses the following properties : 

(1) Any configuration of the assembly can be described by specifying 
the number and geometry of occupiable sites in the assembly and how each 
site is occupied. 

(2) Interactions between pairs of molecules that are not nearest 
neighbours may be neglected, and therefore the interaction energy of the 
systems composing the assembly can for each configuration be expressed 
as the sum of contributions from each pair of immediate neighbours. 

Because of their simplicity a great deal of attention has been devoted 
to the study of regular assemblies. 

Let us first consider a regular assembly consisting of a mixture of two 
types of system which we label 1 and 2. Let V, be the number of systems 
of type 1, N, the number of type 2, the total number of lattice points 
being V (=N,+N,). 

When, the systems are arranged in some configuration on the lattice, 
the energy will be determined by nearest neighbour interactions, and 
let Nj, Nig, Nog denote the number of 1-1, 1-2 and 2-2 pairs in this 
configuration; let €;, €9, €:2. denote the corresponding interaction 
energies. Then the energy of this configuration, is 


N41€41 + Ny9€19 + Nasese. . . . . . . (1) 


However, the numbers N,,, Ny., No. are not independent. Let q be the 
coordination number of the lattice, and let us imagine that each 1 system 
has q arrows attached to it. Then the total number of arrows will be qN,, 

and each 1-1 pair will correspond to two arrows, and each 1-2 pair to one. 


Thus 
and similarly AN sea GN, \ ema arit ae Dp 
Ny+2No,= Ny. J 
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(This is a valid approximation if boundary effects are ignored. Alternatively 
the equations are exact if the lattice is wrapped on a torus. The latter 
device is closely analogous to the use of the cyclic boundary condition, in 
lattice dynamics, and should be asymptotically equivalent to a physical 
lattice for large N.) Hence for given N only two of the five quantities 
Ny, Nz, Ni, Nis, Noo are independent, and we may conveniently choose 
them to be NV, and N,,. The energy (1) of our configuration is thus 


E(N ; Ny, Ny) = BON €n9 + $9.N 1 (€11 — €o9) + Nyelere— $(€11 + €99)]. - (3) 


The final energy in brackets in (3), €1.—4(€,; + €:9), is of basic importance 
in the development, and we shall denote it by «. The energy of any 
configuration with given N,, N,. can therefore be readily determined. 

The fundamental problem in the theory of regular assemblies is the 
determination of the combinatorial factor, g(N ; N,, N,.), which specifies 
for any given lattice the number of possible configurations corresponding 
to given values of N,, Nj. This problem can be formulated very simply, 
but its solution is extremely formidable. In fact only a particular class 
of solutions for two-dimensional lattices are so far available, and even 
these require sophisticated mathematical methods. 

Instead of dealing with g(N ; N,, Nj.) it is more convenient to consider 
the generating function 

ae g(N ; Ny, Nyx) exp (— &/kT) 


= Zt 3 g(N ; N,, Np) wNiz% ° * A . e 5 (4) 
Ni, Nie 
where 


w=exp [—$9(e— €19)/kKT], z=exp(—€/kT), 
29 = OXP (— €y9/2KkT'). 


This generating function will be interpreted physically as a partition 
function, but its detailed character will depend on the nature of the 
component systems. These systems may be atoms or molecules of two 
different types, in which case we have a simple model of a solid solution ; 
Zy will be a type of grand partition function, since we are summing over 
all N,, i.e. all possible concentrations. One of the species may correspond 
to the absence of an atom, in which case we have an assembly in which 
there is a certain degree of free volume, and this enables us to construct 
a crude model of a liquid. Alternatively the systems may consist of 
atomic spins, each spin being capable of two orientations, and we have 
a simple model of a ferromagnetic ; in this case all values of Nj are possible, 
and Zy represents an ordinary partition function, the free energy being 
related to Zy by 


Fy= —kT In Zy = 39QN €:9— kT In A y(p, 2) ne, reel) 


where Ay represents 
SGN as, N49) 22%, 
Ni, Nis 
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We should expect the free energy to be an extensive property of the 
assembly and hence proportional to N for large N; thus on physical 
grounds A, is asymptotically of the form [A(y, z)}” where A does not 
involve V. Similarly we should expect g(NV ; NV, Nj.) to be given in terms 
of the entropy of the assembly by Boltzmann’s relation 


S=kIng(N; N,N»). usa! aa 


Since the entropy is also an extensive property, the asymptotic form of 
g(N;N,,Nys) is f(,y)" where a,=N,/N and y=N,,/N. Quantities 
independent of NV such as 

Z= Z,2A(u, Z), 

F = }4€o2— kT In A(p, 2), } 


will be referred to as the partition function or free energy per system. 


(7) 


2.2. The Ising Model 


It is well known that a non-interacting assembly of atomic spins 
gives rise to paramagnetic properties which vary smoothly and continu- 
ously with temperature. The physical properties of a ferromagnetic 
show discontinuous behaviour at a temperature 7’. known, as the Curie 
temperature. At temperatures below 7’. a ferromagnetic possesses a 
spontaneous magnetization, the magnetization tending to a non-zero 
value J, as the field is reduced to zero. J, is a continuous function of 
temperature which becomes zero at 7'¢; at temperatures above 7’, there is 
no spontaneous magnetization, but the substance possesses a paramagnetic 
susceptibility which becomes infinite at 7c. The specific heat of a 
ferromagnetic has a A-type singularity at the Curie point. 

In, fact the above properties of a ferromagnetic were accounted for in 
a rough way many years ago by the Weiss assumption of an internal 
field. But the origin of this field was not explained, and no attempt was 
made to relate the field to individual atomic interactions, and to solve the 
corresponding statistical problem. 

In 1925 Ising considered an elementary model of a ferromagnetic. 
He assumed that each atom possessed a spin and hence a magnetic moment 
which could orient either parallel or anti-parallel to an external magnetic 
field. There is an interaction between nearest neighbour spins in the 
lattice, parallel spins having an interaction energy —J, and anti-parallel 
spins +J (J>0). Ising hoped that this interaction might be adequate 
to account for the basic properties of a ferromagnetic. In fact he was 
able to solve the corresponding statistical problem in one dimension 
only, and in this case the assembly did not exhibit ferromagnetism. 
Subsequently Peierls (1936) was the first to show rigorously that a 
two- or three-dimensional model would exhibit ferromagnetism. 

It will be seen that the Ising model is an example of a regular assembly, 
although in the presence of an external magnetic field H the interaction 
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of each magnetic moment with the field must also be taken, into account. 
Let ~ be the value of the magnetic moment. Then if the direction parallel 
to the field is denoted by 2 and the opposite direction by 1, the energy of 
an (NV, N,,) configuration, is given by 


(Ny — Ny) mH — (Ny +No)I + Nyy J =(2N,—N)mH—}gNJ+2N J. (8) 


Hence by comparison with (3) and (4) we can write the partition function 
for the Ising model in the form 

Lig ee pie REA (16, Zz), ed ee EY cal 9) 
where now p=exp(—2mH/kT) and z=exp(—2J/kT). The lowest 
energy state corresponds to V,=0, NV,.=0, Ay=1, since J>0. 

In the absence of a magnetic field ~=1. It should be noted that in 
this case there are two equivalent lowest energy states consisting of all 
the spins pointing in direction, 1 or in direction 2 (fig. 1). In the presence 
of a magnetic field this degeneracy is removed, and the lowest energy 
state consists of all spins pointing in direction 2. This point will be further 
considered when, we discuss long-range order in, § 3.6.4. 


Fig. 1 
(ae Us OE Se eae ea Yo te yy OF YY 
NS ie ae a ea ee ee a 
hee eet te St as sae kN ag! 
(a) (0) 
Direction 1 t Direction 2 | 


Ising model of a ferromagnetic. 
In the absence of a magnetic field (a) and (b) are equivalent lowest energy 


states. In the presence of a magnetic field in direction 2 (5) is the only 
lowest energy state. 


The Ising interaction can be simply expressed in terms of the spin 
variables of the individual atoms, assuming this spin to have the value 4. 
Then if s, is the z component of the spin which can take the values + 3, 
the interaction so far discussed between two nearest neighbour atoms 7 
and j is equivalent to 4Js,,s,;. This leads us to a natural generalization 
of the model to an arbitrary spin s for which s, can take (2s+ 1) values. 
The statistical problem now arising is that of a regular assembly of (2s + 1) 
components, although the interactions between, these components bear a 
special relationship to one another. In assessing the effect of a change 
of s, and particularly in allowing s to become infinite (corresponding to 
a continuous range of possible orientations of the spin), it is convenient 
to normalize so that the maximum and minimum values of the interaction 
remain, equal to +J. The generalized Ising interaction for spin s is then 
$48 


a Barer 
2 a 
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In the past few years a good deal of attention has been devoted to the 
problem of anti-ferromagnetism, and if we now consider the possibility 
J <0 we are furnished with an Ising model of an anti-ferromagnetic. 
The formulation, of the problem in (8) and (9) is still valid, but for sufficiently 
small H the lowest energy state corresponds to the maximum possible 
value of N,,. The value of this maximum depends on the nature of the 
lattice; if the lattice can be decomposed into two equivalent sub-lattices 
(e.g. the simple quadratic (S.Q.), simple cubic (S.C.) or body-centred 
cubic (B.C.C.)), the maximum corresponds to a perfect alternation, one 
of the sub-lattices containing spins pointing in direction 1, and the other 
containing spins pointing in direction 2. The maximum value of Nj, is 
then 4qN, and for this state the magnetic field terms in (8) make no 
contribution, to the energy. However, for large H the state in which all 
the spins are aligned in, direction 2, so that VN, =0, N,,=0, will have lower 
energy. The critical field, H., at which the two states have equal energies 
is given by wH,.=—qJ. For fields H< He. there are two equivalent 
lowest energy states, since 1 and 2 can be interchanged (fig. 2); for H > He, 
however, there is only one lowest energy state (fig. 1 (6)). 


Fig. 2 
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(a) (6) 
Ising model of an antiferromagnetic. 


(a) and (b) are equivalent lowest energy states in the presence of a magnetic 
field H < Ho. 


For lattices which cannot be decomposed into two equivalent sub- 
lattices (e.g. the triangular or face-centred cubic (F.C.C.)) the problem 
of the lowest energy state is more complicated, and will be discussed 
further in § 3.5.4. 


2.3. Regular Solutions 


A rich variety of physical phenomena is associated with binary mixtures, 
many of them, such as the separation into two phases, corresponding to 
a discontinuous behaviour of the thermodynamic functions. In considering 
the properties of binary solutions we shall again use the simplest possible 
model capable of accounting for the basic facts. These are that under 
certain, conditions two substances can mix freely at all concentrations, 
whereas under other conditions, below a certain critical temperature 7'o, 
separation into two phases occurs at definite concentrations. Thus there 
is a solubility curve for any one substance in another which is a continuous 
function of temperature. 
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We shall assume that our binary mixture satisfies the conditions laid 
down in § 2.1 for regular assemblies, and we shall refer to solutions of this 
type as ‘regular solutions’. This term was first introduced by Hildebrand 
(1929) in his pioneering work on the properties of non-ideal solutions. 
Hildebrand defined a regular solution as one in which the mixing of the 
constituents was completely random although the energies of interaction 
between different molecules were different. This concept proved very 
fruitful in a semi-empirical analysis of the properties of solutions. 
However, from the theoretical point of view, this can only be regarded 
as a first approximation, since the positions of molecules in a solution 
must be correlated with the interaction energies. For this reason Fowler 
and Guggenheim introduced the alternative definition quoted in § 2.1 
which provides a simple model from which the properties of solutions 
can be calculated in a fundamental way. We shall always use the term 
in this sense, as it has become widely adopted for theoretical work (see 
also Hildebrand 1951, 1953, Guggenheim 1953). 

In considering the partition function for a regular solution of given 
concentration the numbers N,, VV, of systems of each kind are fixed and 
we do not sum over these variables. Using the energy interactions of 
§ 2.1 we may write for the partition function for a regular solution 

Zy§ = ees > g(N ; N,, N,,)2"™, oi Mabe” 9) ce (10) 
where 2,=exp(—«,/2kT), z,=exp (— €9./2kT), z=exp(—<«/kT). If we 
denote, as before, the concentration of component 1 (=N,/N) by a, 
then, as V becomes large we can define a partition function per particle, 
independent of NV, given by 


ZB = 210%), h( cx, 2). jae Pee (ET) 


It will be seen that the sum over J,, in (10), which is asymptotically 
equivalent to [s(a,,z)]}", is equal to the coefficient of 11 in Ay(p, 2) 
from (4). Hence we may write this as a contour integral 


al hd 
= Ay(p, 2) dy [pre 
74} Co 


where C is a circle surrounding the origin. For large N the integrand 
has a saddle point on the real axis if there is a value of » satisfying 
(d/dy) (NInA—N,Inp)=9 
or 
(wd/dy) (In A) = oy. Pee bind al 2) 


Following the usual argument, by choosing our circle to pass through 
the saddle point, the major contribution to the integral comes from the 
neighbourhood of the saddle point, and we may write asymptotically 


In f(x, 2) =In A(u, z)-—% Ing, Mac (13) 
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where ; is now a function of a, and z given by (12). Hence if we know 
A(p,2) as a function of the two variables Hy & and we wish to determine 
the regular solution function x(a, 2) for a given concentration «,, we must 
first determine the curve in, the (u,z) plane given by (12), and we then 
use (13) to obtain y(a,,2). Thus the regular solution function is completely 
determined by a knowledge of A(x, 2). 

The result (13) could equally well have been deduced by replacing the 
sums in (4) and (10) by their maximum terms. A more fundamental way 
of understanding the significance of the result is to derive it from the 
grand partition function (Rushbrooke 1949, Ch.17). For a regular 
assembly consisting of two components the grand partition function is 
given by 

B(T, Ay As) = > Agen) Agee Nal a2 


Ni, Ns, Nv 


and hence from (4), this can be put in the form 


E(T, Ay, A2) = DY (Avz24)*An(u, 2), er 0), 

N 
where pw is now equal to A,2,7/A,z.2.. Since we are concerned with large 
values of NV we can replace Ay by its asymptotic form A”, and we find that 
B= 1/(V—Agegth) eS a eee 


Using the standard formulae for the G.P.F. 


N, =A, 0/0, NE =D,z42 A/Au/(1 —Agza2A), an 


Ny=z 0/0Aq In E = (Agzo27A — A424 OA/Opn)/(1 —Agzo%A). 
Therefore, the concentration « of component 1 which is equal to N,/N 
or N,/(N,+N,), is given by 


het OA 8 
AgzotA On “ne Op 


My se 


(hiA);, “. 8 \se So ceee 


and this is in agreement with (12). 

It will be seen from (15) that the function A,(u,z) which arose in the 
Ising model can be regarded as a ‘restricted’ grand partition function 
in which the ordinary partition function 1s summed over all possible 
concentrations, but the total number of systems remains fixed and equal 
to NV. The thermodynamic properties of regular solutions can readily 
be determined from the G.P.F. (16). 

We must now examine the lowest energy state of a regular solution of 
given, concentration, and we see from (3) that this depends on, the sign of e. 
If «>0 the lowest energy state corresponds to minimizing N,,, and this 
can, be achieved by dividing the assembly into two separate sections of 1 
and 2 atoms with the minimum region of contact. There may be a number 
of independent equivalent ways of choosing this minimum region, but 
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they will not give rise to any entropy. We may say that the lowest state 
of the assembly consists of two separate homogeneous phases. ‘This 
sign of € corresponds to the Ising model of a ferromagnetic. 

If however « <0 the lowest energy state of the assembly corresponds to 
maximizing N,,, which can be achieved by surrounding each 1 atom by 
as many 2 atoms as possible, and vice versa. This results physically in 
the formation of a superlattice, and hence the negative sign of « is 
associated with order—disorder transitions in alloys; it also corresponds 
to the Ising model of an antiferromagnetic. The lowest energy state is 
degenerate, and except for a particular concentration, for a given lattice 
gives rise to a finite value of the entropy. The determination of this 
entropy, and of possible lowest energy states is quite complicated, and is 
related to the anti-ferromagnetic in a critical field. 

The problem of regular solutions with more than 2 components is much 
more involved. Formally a relationship with a corresponding Ising model 
can, be established and results such as (12), (13), (15) and (16) can readily 
be generalized. However, the interactions between atoms in various 
possible spin states must be arbitrary for complete generality, and the 
complications arise from the different possibilities of magnitude and sign 
of these interactions. 


2.4. Surface Problems 


Another group of phenomena in which interactions between elementary 
systems can play an important role are the adsorption, of gases on, solids, 
and the growth properties of crystal surfaces. Langmuir proposed a 
model of an adsorbing surface as consisting of a fixed number of sites on 
which adsorption, can, take place (see Miller 1949, Ch. 2). Each gas 
molecule that is adsorbed occupies one site, and the sites are arranged 
in a regular two-dimensional array. If we assume that forces between 
molecules occupying sites which are not nearest neighbours can be 
neglected, our model is another example of a regular assembly. The 
forces between nearest neighbour molecules may be attractive or repulsive, 
and we wish to determine how the heat of adsorption and the fraction, of 
surface covered depend on the pressure and temperature of the gas 
adsorbed. We must assume that thermal equilibrium is attained in the 
adsorbed film, and we are therefore concerned with what is usually called 
a mobile film. 

The systems 1 and 2 of §2.1 now correspond to occupied and vacant 
sites respectively. Hence we now have ¢,.= €92=0 in (3), and «= — 3¢). 
However, each adsorbed atom interacts with the surface as well as with 
its nearest neighbours; if we denote the energy required to evaporate 
one adsorbed atom from an, otherwise bare surface by 7, the energy of an 
(N ; Ny, N42) configuration is , 


—Nyyt+N yp. A e ° ° . e 4 (19) 
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Hence from (4) the partition function for an adsorbed monolayer of N sites 
is Ay(,2), Where now »=exp(n/kT). The fraction of occupied sites is 
given, as in $2.3, by (12). The conditions of equilibrium between the 
layer and the adsorbing gas can be determined in the usual manner by 
equating the thermodynamic potentials in the two phases. 

Similar considerations apply to the equilibrium structure of a crystal 
surface, a problem which has been discussed in detail by Burton et al. 
(1951). The ‘adsorbed’ layer now consists of the same molecules as the 
crystal, and we are interested in properties such as the rate of growth 
and surface roughness as functions of temperature. If we assume that 
only a single layer of molecules can be added the treatment is the same 
as for an adsorbed monolayer except that the atoms of the monolayer are 
the same as for the solid. If we allow for the possibility of more than one 
molecule becoming attached to a given site, and also allow molecules to 
become detached leaving vacancies, the resulting problem is somewhat 
similar to the generalized Ising model with arbitrary spin. However, 
the interaction between different levels is different, and is proportional 
to |s,;—s,| where s; and s, take positive or negative integral values or zero. 


2.5. Conclusion 


It will be seen that the properties of a regular assembly with two 
components are completely determined by a knowledge of a function, of 
two variables A(,z) which is the partition function, per spin of the Ising 
model in the presence of a magnetic field. There will be one particular 
function, A(,z) corresponding to each type of crystal lattice. Once the 
nature and position of the singularities of the Ising model are known the 
physical results can be interpreted for any of the other problems discussed 
above. Our first task, therefore, is to obtain as much information 
as possible about A(u,z) for different types of lattice. We need consider 
only values of «<1 since values of «21 merely correspond to an 
interchange of the labels 1 and 2. However, the values of z can range 
from 0 to oo, and there is a fundamental difference between the case 
0 <z<1 (e>0) which corresponds to a ferromagnetic, and 1 <z< 0 (e<0) 
which corresponds to an antiferromagnetic. 

The exact evaluation of A(u,z) for all « and z, has so far been possible 
only for a one-dimensional model. For a number of lattices in two 
dimensions A (y,2) has been evaluated exactly when = 1, and some other 
properties of the model have been determined exactly. In three dimensions 
no exact information is at present available. However, several general 
theorems have been established which enable us to say a good deal about 
the position of the singularities. In the absence of exact solutions two 
approaches have been used. Expansions have been developed as power 
series in yw, z, and 1/z, 1/7’ (near z= 1); and various approximate evaluations 
of the function have been undertaken, which correspond to ignoring 
certain features of the exact solution. We shall now proceed to discuss 
these points in further detail. 
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Notation for § 2 


N total number of units in assembly. 
q coordination, number of lattice. 
N,,N, number of units of types 1, 2 in assembly. 
a4 =N,/N, a=N,/N. 
Ny, N2,Ny, number of nearest neighbour pairs of various types. 
E11 €99, 2 energies of interaction of nearest neighbour pairs of 
various types. 
€ = €19— 3(€11 + € 9). 
% =exp [—(<,/2k7)], 
%=eXDp | — (€g9/2kT)]. 
z=exp [—(e/k7’)]. 
T =absolute temperature. 
k= Boltzmann’s constant. 
7T'.= critical temperature. 
g(N ; Ny, Nj.) =combinatorial factor, i.e. number of configurations with 
specified values of V,, N49. 
Zy,F'y,Ey,Syty partition function, free energy, internal energy, entropy 
and magnetization, of assembly of V units. 
Z,F,H,S,I corresponding functions per unit of the assembly. 
Ay(p, 2) = > yN; Ny, Ny_)p%12%e ; he [Ay]! = A(p, 2). 


I> N12 


Ising Model. J=interaction energy of parallel spins. 
H =external magnetic field (taken, in direction 2). 
m=magnetic moment of individual spin. 
z=exp[—(2J/kT)], w=exp[—(2aH/kT)]. 
s=total spin; s,, s,, 8, components of spin. 
H,=critical field for an anti-ferromagnetic. 
Zy', Z| partition function for Ising model. 
Z,§, ZS partition function for solution. 
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§ 3. Exact Mretuops AnD SOLUTIONS 
3.1. One-dimensional Assemblies 


For a one-dimensional assembly of systems each having ¢ possible 
energy levels a general treatment is possible which reduces the evaluation 
of the partition function to the determination of the largest eigenvalue of 
a matrix. Let U(x;,x;,,) denote the energy between any two systems x; 
and 2;,, each x; taking ¢ possible values. We shall assume that the systems 
are connected in a ring so that the mth system is connected to the first ; 
this considerably simplifies the algebra and is analogous to the cyclic 
boundary condition in lattice dynamics. Then the partition function 
for m systems is given by 


Zin= >, XP i = [U (ay, %2) + U (ag, %3) +... U(%__1, Um) + U (%m; “lf 


t 
SD Vig Viewed aa es Ss Oe ee 
i, jud=1 

where V,,=exp(—U,,/kT’) corresponds to the interaction between 
neighbouring systems in states 7 and j respectively. The elements JV;; 
thus form a txt matrix V and the continued product (1) will be readily 
identified as the trace of VW". Thus if A,,...A,;f are the latent roots of V, 
Zp =m tA +... FAM. n'y, ty eel eee 
For sufficiently large m only the largest root A, is significant, and we may 
identify Z,, with A,”, and the partition function per system with j. 
Even, if the largest root is p-fold degenerate this result still holds since 
pm tends to 1 as moo. It should be noted that if the interaction V;,; 
is symmetric in 7 and j, then the matrix V is symmetric; also (V”),,, the 
ith diagonal element of V” is proportional to the probability that a 

particular system is in state 7. 

The above method of treating a class of one-dimensional assemblies was 
suggested independently by Kramers and Wannier (1941), Lassettre and 
Howe (1941), and Montroll (1941). It can be applied to systems with a 
continuous range of energy levels, but the eigenvalues in (2) must now be 
replaced by the eigenvalues of an integral equation 


| expt- Uly,2)/kTJaly)dy=Aa(z). . . . . (3) 


In this case the number of eigenvalues will no longer be finite and additional 
care with the pure mathematics is needed to establish the conclusions. 


3.2. Application to the Ising Model 
3.2.1. The simple Ising chain 


For the simple one-dimensional Ising model in the absence of a magnetic 
field the application is immediate. There are two states for each system 
and V is a 2x2 matrix; for a generalized r-state model V would be an 


+ These should not be confused with A,, A, of § 2.3. 
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rxr matrix. However, in the presence of a magnetic field, the interaction 
with the field must also be taken into account, and to preserve the 
symmetry of the treatment it is sometimes convenient to attribute half 
of the interaction of each of the two neighbouring systems x; and x;,, to 
U(x; %;,,). We thus write for the generalized Ising model of spin s 


mH 
gi 5082 it1~ —o— (82,4 +Se,441)) Homes tey we) 


U (4,504) = — 


12 t-- 1 


and in particular for the two-state model, using the notation of § 2.2: 


V == —1/2,—1/2 1/2 
: 1/2 :. 1/2 (5) 
z 2 wi 
and 
Awz)\=HO+m)+-VI—pP +422. 8... (8) 


The two eigenvalues A,, A, correspond to positive and negative signs of 
the root. 

Despite the simplicity of this result a number of interesting and useful 
conclusions can be drawn from it. In the absence of a magnetic field (1 = 1) 
the partition function reduces to z—¥?+z!2 (see (9), §2). There is no 
singularity in the partition function, which is equivalent to that of an 
assembly of non-interacting systems each having two independent energy 
levels. The specific heat is equal to mk(J/kT')*sech?(J/k7’), and this is 
drawn, as a function, of temperature in fig. 1. Although the specific heat 
has a maximum at k7’=0-833 J this maximum is not associated with any 
peculiar behaviour of the magnetic properties of the chain. 

In fact the magnetization of the chain per spin is given by 

L=kT 0/0H [In Z] = —2mp0/ou.[—4lnp+m A] 
= m[1— 2uA,/A]. Seecog okt) 
Now, 

—2(1—p) + 42? (8) 
4y/[(1—p)? + 4y2"]’ 
and hence as H-0, and p>1, A,>}(1+z2), and I-+0 for all z. Thus the 
magnetization tends to zero as the field tends to zero for all temperatures, 
and there is no spontaneous magnetization. However, if we proceed to 
determine the susceptibility we do find evidence of peculiar behaviour at 
T'=0. The susceptibility is given by 


A,=4+ 


4n? 
x= O1/0H = —2nm|kT OL/Ou= : - [wA,/A+p2(A,,/A—A,2/A?)], 
(9) 
and as p1, this is equal to 
(4o2/kT')[4 BA alk =m|kTexp(2J/kT). «© « . (10) 


Thus when J>0 (i.e. for a ferromagnetic) the susceptibility becomes 
infinite at 7’=0. If we now examine the expression (6) we see that a power 
series expansion, in ». and z is possible for A(u, 2) about any point for which 


P.M.S. % 
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u#1, 240. No such expansion is possible, however, about »=1, z=0, 
which corresponds to zero temperature and zero magnetic field. For this 
reason, it is often, convenient to describe a linear chain as a ferromagnetic 
whose Curie point has gone to zero. 

We can, also learn, something of the behaviour of an antiferromagnetic 
from the solution, for the linear chain. In the absence of a magnetic field 
the partition, function, is the same as for a ferromagnetic. However, the 
magnetic behaviour is very different. The susceptibility in zero field is 
given by (10) but J is now negative; hence it tends to zero as 7'>0 and 
as 7->0oo and has a maximum at a temperature 2J/k. The behaviour at 
T' =0 is no longer peculiar in zero field but in a critical field Hp=2J/~ 


Fig. 1 


5 
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Specific heat of a linear Ising chain. The maximum does not correspond to 
anomalous thermodynamic behaviour. 


(see end of §2.2). Let us consider the lowest energy antiferromagnetic 
configuration, which consists of alternating spins on each of two sub- 
lattices pointing respectively in the direction of, or in the opposite direction 
to, the field. In the critical field any spin of the latter sub-lattice can be 
overturned without change of energy. Hence the lowest energy state has 
degeneracy at least 20", But in fact there is further degeneracy in the 
assembly, since if any two neighbouring spins of the second sub-lattice 
have overturned to point in the direction of the field, the spin of the first 
sub-lattice which lies between them can also be overturned with no 
change in energy (fig. 2). We can calculate the actual degeneracy in the 
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critical field from (6), remembering that when »/] = —2., p2z2=1. We 
have, for the entropy, 


S=—0F/0T=kinA+kT 0/dT (In A), soe e (lal) 
and, when ..z*=1 the second term gives zero as T->0. Hence 


Se>kIn (- ) =kIn (1-618). 


Fig. 2 
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Lowest energy states of an antiferromagnetic chain in a critical field. (a) Any 
of the underlined spins can overturn without change of energy. (b) 
if two neighbouring underlined spins in (a) have overturned the inter- 
mediate ringed spin can overturn without change of energy. 


Fig. 3 
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Entropy of an antiferromagnetic chain as a function of magnetic field. The 
initial field breaks up the alternating order and increases the entropy. 
When the field increases sufficiently ferromagnetic ordering takes place, 
and the entropy decreases. 


N2 
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At 7 =0 if H is slightly less than, or slightly greater than H, the entropy 
is zero again. Thus the entropy as a function of magnetic field consists 
of an isolated line at H=H,. At temperatures different from zero this 
line spreads out into a curve which becomes broader as the temperature 
rises, but the maximum is still recognizable at quite elevated temperatures 
(fig. 3). This corresponds to the physical fact that a magnetic field breaks 
up the alternating order of the spins, and hence increases the disorder of 
the chain; only when the field increases sufficiently to line up the spins 
does the order again increase. As might be expected the susceptibility 
in the critical field also has anomalous behaviour at low temperatures ; 
it is easy to show from (9) that y>/5m?/25kT. 

The low temperature behaviour discussed in the previous paragraph 
is clearly in disagreement with the third law of thermodynamics; such a 
system giving rise to a finite entropy at the absolute zero for a particular 
magnetic field could be used to reach absolute zero (Simon 1956). From 
this we may conclude that the Ising interaction is not a rigorously correct 
form of physical interaction, and for a real physical system additional 
interaction, terms would enter to split the degeneracy in a critical field. 

It would be of interest to investigate the corresponding behaviour of. 
a linear chain with a Heisenberg interaction, Jo;.o; between neighbouring 
spins (where the o; are the Pauli spin operators) since this interaction 
has a sounder physical basis. Although a formal solution exists to this 
problem (Bethe 1931, Hulthen 1938) it has not so far been possible to 
use this solution to calculate physical properties such as those discussed 
above. The most general interaction between spins which takes account 
of magnetic anisotropy is of the form Ao,0,;+ Bo,o,;+Co,0,;, and the 
Ising interaction is thus an approximation corresponding to extreme 
magnetic anisotropy (A and B negligible compared with C). For certain 
substances (e.g. FeF,) this can be regarded as a reasonable approximation ; 
however, at sufficiently low temperatures the terms A and B can no longer 
be neglected. 

Before leaving the one-dimensional solution we should draw attention 
to an interesting difference between the ferromagnetic and antiferro- 
magnetic solutions when, m is odd. The exact solution for the partition 
function in the absence of a magnetic field is, by (2) and (6) 


Ly = (2-H? + gu2ym 4 (g-W2—gi2ym, | (12) 


When m is even this expression is invariant when J is replaced by —JJ. 
However, when m is odd the second factor is multiplied by —1, and this 
corresponds to the physical fact that one spin must be a misfit in the 
ordered antiferromagnetic chain. Thus the lowest energy state for the 
antiferromagnetic case has energy (m—2)J because of the energy of the 
misfit spin, and is 2m-fold degenerate because of the position which this 
misfit can occupy. However, at any temperature different from zero the 
second term in (12) is negligible compared with the first for large m, 
and the ferromagnetic and antiferromagnetic solutions are identical. 
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Thus the misfit spin contributes to the partition, function only at the 


absolute zero. We shall pursue this point further for the two-dimensional 
model in § 3.5.1. 


3.2.2. The two- and three-dimensional models 


The same method can be applied to the two- and three-dimensional 
Ising models, building up the lattice by adding one row or one layer at a 
time. Thus for a two-dimensional lattice consisting of n rows as shown 
in fig. 4, each system of §3.1 consists of a row of n spins with 2” possible 
energy states. Thus the matrix V(= V,,) is of order 2” (r” for the generalized 
model). The treatment of § 3.1 will apply if the systems are joined to form 
an anchor ring. Similarly for a three-dimensional lattice consisting of 
p xn layers joined together to form a torus the matrix V is of order 2”. 
It should be emphasised that the matrix V(=V, ,,) can be chosen to be 
symmetric, and all its elements are positive. 


Fig. 4 
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Each added unit consists of a row of ” spins. 


However, it is quite a complicated matter to write down the appropriate 
matrices V,, V,,, for particular values of p and n, and these matrices 
can be decomposed into product matrices of simple characteristic structure. 
This can most readily be demonstrated by the use of the ‘screw trans- 
formation,’ introduced by Kramers and Wannier (1941). In the two- 
dimensional lattice consisting of n rows we now build up successive rows 
by adding one unit at a time instead of a whole row, as shown in fig. 5. 
The treatment of §3.1 can be modified to deal with this, although the 
symmetry of the matrix must now be sacrificed; V,, can be written in, the 
form A,,”-1B,,” where A,,, B,, have simple characteristic structure. 

The B,, matrix corresponds to an end effect, and if successive columns 
are linked together to form a screw, the nth elements of any row being 
linked to the first element of the succeeding row, only the matrix A, 
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remains. It seems very likely that as n+ o the effect of the screw linkage 
will become negligible, and the partition function for the two-dimensional 
quadratic lattice will be given by the largest eigenvalue of A,. Ananalogous 
transformation can also be carried out for the three-dimensional model 
(for further details see Domb 1949 a). 

Despite the apparent simplicity of structure of these matrices little 
progress has been made in reducing them or determining their eigenvalues 
except for two-dimensional lattices in the absence of a magnetic field 
where other solutions are also available (Newell 1950). They were instru- 
mental in leading Kramers and Wannier (1941) to the exact location of 
the Curie point of the simple quadratic lattice, and they have been used 
to derive power series expansions for the partition function (Ashkin and 
Lamb 1943, Domb 1949b). However, these can usually be obtained more 
easily by combinatorial methods which will be described in § 5. 


Fig. 5 
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Kach added unit consists of a single spin. 


3.2.3. Conditions for the existence of singularities 


The matrices V introduced in the previous section are finite and 
symmetric; their elements are simple analytic functions of ~ and z 
(u, #0) which are bounded and have bounded derivatives. To determine 
the derivatives of the largest eigenvalue in the neighbourhood of a given 
point we can use standard perturbation theory as developed for example 
m quantum mechanics (e.g. Mott and Sneddon §18.1). At any given 
point the matrix has eigenvalues ),,...A, and corresponding eigenvectors 
E1,---§, and even if the eigenvalues are not all distinct, because of the 
symmetry of the matrix the eigenvectors can be chosen to be independent 
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and orthogonal. The changes in A, and €, as a result of a change 8V in V 
are given by 


dA, = er eV En + second-order 
, had bey. i 
0b1= & 51 OVE; /(41—A,), | terms in dA, 68. ae 


Thus if the largest eigenvalue is non-degenerate both $A, and S€, will be 
of the first order, from (13), and the limits to which 8), 1/Ou, 5A,/dz tend 
will be independent of the paths of du, 5z. Hence A, will be an analytic 
function, without singularities. 

An alternative method of deriving this result is to consider directly 
the characteristic equation for V. It is easy to see that this is of the form 


Oa ire -> Nobi( WL, 2 Ngee PS ONT Dee OE 
where ¢,(u, 2) are simple analytic finctions of » and z (u,z40). Hence 

0A, /0z= — (Of/dz)/(Af/OA),_), ; tat hla (LO) 
A, will only be singular if 0f/0A=0, ie. if A=A, is a repeated root of 


f(A; pt, 2) =0. 

Now there is a theorem of Perron (1907) and Frobenius (1908) which states 
that finite matrices with positive elements have largest eigenvalues which 
are simple. From this we draw the conclusion, that any lattice which is 
infinite in one dimension only will have no singularities. The argument 
can, readily be generalized to any forces of finite range. A singularity can 
only arise for such forces if the lattice is infinite in at least two dimensions. 

The proof given by Frobeniusf is involved and depends on the use of 
mathematical induction. In our applications only symmetric matrices 
need be considered, and a simpler proof can be constructed. Since the 
theorem is of some significance in the theory of cooperative phenomena 
we give such a proof in Appendix I. 


3.3. Order and Long-range Order 

Let <a,) and <a.) represent the mean fractions of spins pointing in 
directions 1 and 2 respectively at any given temperature. Consider any 
two sites in the lattice which we denote by a, 0; it is of interest to know 
whether the spins on, these sites are correlated, and we wish to find a 
method of defining the correlation precisely. Let (a4) (a9) (oy) 
(o>, represent the mean fractions of pairs of spins on, sites equivalent 
to a and b being occupied at any given temperature by 1-1, 2-2, 1-2, 2-1 
respectively. Then it is clear that the following relations hold 


{041% ) + Cay2™ ) = {x ), 
{Ola ) + Cotgn7?) = Cola), Pie ee UO) 
Cayo") = Cola”), Cat + Cog) = 1. 
Hence if (a,) and (a ,) are given (a417) Cayo”) Cog”) depend on only 


+ Alternative proofs have been constructed recently—see Gantmacher (1959), 
Bellman (1960). 
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one parameter, and this can be used to define the order between the sites 
a and b. If there were no correlation between the spins on a and 6b we 
should have 


Coty?) = (ax )®, Cayo”) = Cay) (He) Cotgn)) = Cay). (17) 
Thus in general we shall write 


(0447?) = (ay )>(1 + Pap), 
and hence from (16) 
Cory) = Cary D[ (2) — C041 DP a]; <i ni g- Le 
(ag0"” ) = Calg) + Coy Pa 
where for convenience (a; ) is taken <}. Ifc¥,,=0 there is no correlation 
between, the spins, and the deviation of c”,, from zero provides a measure 
of the correlation. 

We now enquire what happens to o¥%,,, as the distance between a and 6b 
becomes very large. If c%,, tends to a limit c¥ ,, different from zero we 
say that long-range order is present in the lattice. To secure agreement 
with the usual definition of long-range order we must write 0% ,, = R®, and 
R is a measure of the long-range order. The advantage of c/ ,, is that it 
can, be both positive and negative, and hence for an ordered antiferro- 
magnetic array its limiting behaviour consists of an alternation in sign 
as the number of steps between the two sites changes from odd to even. 
If we are to use RF, the definition (18) is valid for an antiferromagnetic 
array only for sites separated by an even number of steps, and a re-definition 
is required for an odd number of steps. 

The above definitions of order and long-range order are immediately 
adaptable to solutions and superlattices containing two components of 
fixed concentrations. If more components are present a number of 
additional parameters are required to specify the order of the lattice. 
This definition, has the advantage that if does not refer to any particular 
basic arrangement of spins on the lattice; hence the concept of ‘right’ 
and ‘wrong’ atoms is not required, and equivalent ordered states are 
automatically taken, into account. 

For the Ising model of spin 4 it is convenient to introduce variables o; 
for each lattice point 7 where the z-component of spin, s,;=40;; o; takes 
the values +1 according as the ith spin is pointing in direction 1 or 2. 
The correlations can then be expressed in terms of the averages <c,), 
{o,0,); it is easy to show that 


(oye Ge etey: 
\ (19) 


(F4Fp) =[ Cot) — Calg) P+ 4404 POF ay 


These formulae can be generalized to provide a definition of order for 
the Ising model with arbitrary spin s in terms of (s,,) and (s,,8,,)3; and 
for the Heisenberg model for which the s,, for individual sites are no 
longer constants of the motion, and one cannot therefore refer to a particular 
spin, as pointing in direction 1 or 2. 


Cooperative Phenomena in Orystals 173 


The concept of order is particularly useful for understanding the 
behaviour of cooperative assemblies. One of the most important physical 
conclusions of our investigation will be that the propagation, of order to 
large distances does not require long-range forces. Even crystals with 
nearest neighbour forces can. manifest long-range order which is then 
destroyed by thermal agitation at certain critical temperatures. 


3.3.1. Condition for long-range order 


We have seen in §3.1 that the partition function of a linear chain, of a 
large number m of systems each having ¢ possible energy levels, is given 
by the largest eigenvalue of a matrix V; and that the probability distri- 
bution of any particular system among the levels at a given, temperature 
is determined by the diagonal elements of VW”. In § 3.2.2 we have shown 
how this can be adapted to the two- and three-dimensional Ising models ; 
for an n-rowed simple quadratic lattice t= 2”, and for a p x n simple cubic 
lattice ¢=2?". We shall now investigate the eigenvectors and eigenvalues 
of the matrix V, and derive their physical interpretation. 

V can be taken as symmetric; therefore its eigenvalues are all real, 
and its left and right-hand eigenvectors are the same. Let us denote the 
eigenvalues by Aj, As...A,, and the normalized eigenvectors by &,...&,; 
the latter can always be chosen to be orthogonal even if there is degeneracy 
in the eigenvalues. Let a denote the vector with a 1 in the ath position 


and zero elsewhere. Then 
t 


a= > a&; Meee ae tee pms eb yah (20) 


t=1 
where «; is the ath component of §;. Hence the ath diagonal element of 
V™, which is equal to a’ Va, is 
t t t 
( >) a6; ( > Ama) => > a 2A,™, ° ° 5 . (21) 
j= i=1 i=1 
and the probability of a given system being in state « is equal to 


t t 
DS ohm | » Ages ° . . ° ° 3 A (22) 
4=1 i=1 


When m is large only the largest eigenvalue need be considered. If this 
is simple, the probability reduces to a,?, i.e. the squares of the components 
of the eigenvector &, are equal to the probabilities of the systems being 
in the state corresponding to this component. If it is degenerate we must 
sum «,2 over all the independent eigenvectors corresponding to the largest 
eigenvalue. For the two- and three-dimensional Ising models, we may 
be interested in a particular spin of a given layer, and we must then sum 
over all states of the layer consistent with the spin having this value. 
We now consider the probability of two systems in layers / and h being 
simultaneously in states « and B respectively. This probability is propor - 


tional to 
t 


> VuaiViy eee Vie a eee Vee ° . e e ° (23) 


Oe yaoi ci 
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and this is equal to the «8 component of VW’ x the Bx components of 
v"-h. Following the same procedure as in (21) we find that «8 component 
of V" is equal to 

t 


Dah Ape we sae 


i=1 


Hence the required probability is 


t t t 
( S aPA" )( s asf" / ¥ jt Pee eee 
i=1 j=l i=1 

When h, m are large the behaviour is again dominated by the largest 
eigenvalue, but we must consider separately the cases when it is simple or 
degenerate. If the largest eigenvalue A, is simple (25) reduces to «,°f,?, 
which is the product of the probabilities for the two layers separately ; 
there is thus no correlation between these layers, and long-range order is 
absent. If on the other hand the largest eigenvalue is degenerate, let 
us suppose, for example, that there are two independent eigenvectors 
which we denote by §, and §_. (25) now reduces to (~,8,+«_$_)? whereas 
the probabilities for the separate layers are a,?+a_”, B,2+8%. We 
thus see that there is a correlation between the layers, and long-range order 
is present. If we wish to evaluate this long-range order for the Ising 
model using the formulae of the previous section for correlations between 
individual spins, we must sum over all configurations of the layers consistent 
with the spins having particular values. We may thus draw the general 
conclusion that the condition for long-range order to be present in an 
assembly is that the largest eigenvalue of the corresponding matrix is 
degenerate ; the long-range order is determined by the largest eigenvectors. 

If we wish to know the value of the order at any finite distance we must 
study the remaining eigenvalues and eigenvectors; for large m all terms 
other than those corresponding to A, in the second bracket of (25) can 
be neglected. The formula for the joint probability of layers 1 and h 
then, reduces to 


auf | Bit ¥ Os)" | . Aaa ee 


if A,, is simple, or 
A. +oP)| oP taf + ¥ Oo]... (7) 


if A, is degenerate. In the first case the ratio (A,/A,) provides a rough 
estimate of the distance to which order extends. We may thus summarize 
by saying that the remaining eigenvalues and eigenvectors of V determine 
the © propagation of order’ in the lattice. This concept was first intro- 
duced by Zernicke (1938), and we have here largely followed the develop- 
ment of Ashkin and Lamb (1943). 
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3.4. Transformation Properties of the Ising Model 
3.4.1. General transformations 


We have seen from § 2 (5) that the partition function for the Ising model 
can be expressed as a sum of terms each of which represents a configuration 
with a given number, N,, of spins pointing in direction 1, and a given 
number, N,,. of 1-2 nearest neighbour links. There is complete symmetry 
between directions 1 and 2; thus if we replace H by —H the partition 
function must be unchanged. Hence from § 2 (9) 


peOPY An (w, 2) = wt? Ay (u-4, 2) 
or Ee te NO) 
pSPA(u, 2) =ptPA(u-, 2). 


If we start with the lowest energy state at the absolute zero, we can classify 
the configurations in terms of increasing deviation from the lowest energy 
state as more spins overturn, and more links are broken. We shall now 
confine attention to lattices such as the S.Q., S.C. and B.C.C. which can 
be divided into two equivalent sub-lattices. In the absence of a magnetic 
field, there is an exact correspondence in energy between the configurations 
of a ferromagnetic (J>0) and those of an antiferromagnetic (J <0). 
This means that the partition function in the absence of a magnetic field 
for these lattices must be invariant under the transformation J—>-— J, 
and hence from § 2 (9): 


gee lee) =e 4A 27), \ 
gaat (lee) = ct A(1 33"). 
Using the o; introduced in the previous section the partition function 
of the Ising model for any lattice can, be written in the form 


Gees {(I@)| Seves+ Fn sath Ses ai], . ey 


where the sum must be taken over 2” possible values of the o; for all the 
N lattice points, and i, 7 are nearest neighbours; we can thus consider 
the first term in, (30) as a sum over all nearest neighbour links in the lattice. 
Clearly the o variables commute and (30) can be written, as a product 


(29) 


Zy' = [[ exp K(e,0; Tex (Fr :) = JV ALN  e eeel Y) 
4,9 
The o,0; eutiate the relation 
(o CO; a)" = (a,0;)4 one = 1, (o,0;) => (o;0;)°= (o;0;)? = eee, 
and hence we can, write 
exp (Ko,o;)=coshK+o,;smhK. . . . . (82) 


We can expand the first product in (31) as follows (van der Waerden 1941) 
TI (cosh K + 0,0; sinh K)= elite +o,0, tanh K] 
ij 
=cosh”K S [1+ (tanh K)(o;0;) 


Weak 


+ (tanh? K)(o;0;)(o,0)) +... J Cane 5. (89) 
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where VW’ is the number of nearest neighbour links in the lattice. This 
can now be interpreted topologically. The first term 0,0; corresponds to 
any one of the ./ nearest neighbour links in the lattice; the second term 
corresponds to any non-identical pair of nearest neighbour links and there 
will be. (MW —1)2! such pairs; similarly for higher terms (see fig. 6). 

The series development (33) can, be used as a basis for a high temperature 
expansion of the partition function in a magnetic field, and will be discussed 
further in § 3.6.2. It should be noted that there is no approximation in 
(33), and provided the sum is carried on to the th term the result is 
valid for any aggregate of spins. 


Fig. 6 
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Typical contributions to the high temperature expansion. 
(a) First term A contributions ; 
(b) Second term AW (4% —1)/2! contributions ; 
(c) Third term W(AW —1)(W —2)/3! contributions. 


We now particularize to the case H = 0, and this introduces a substantial 
simplification, since the second product in (31) becomes unity. Also any 
configuration with an odd number of lines meeting at a vertex will have 
an, odd o; left in the summation, and will give zero contribution. Hence 
the only non-zero contributions arise from closed polygons each vertex 
of which is the meeting point of 2, 4, 6... nearest neighbour links. Typical 
polygons are shown in fig. 7. 

Any such polygon gives a contribution 2 for each vertex when, we sum 
over the appropriate o;; the o, at all remaining vertices likewise give a 
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factor 2, and we obtain a total factor of 2. Thus we can, now write the 
partition function, in the absence of a magnetic field in the form 


MW 
Zy! = 2*(cosh Ky" | 1+ > P(r) tank | Sterna) 
r=1 


where p(r) is the number of independent polygons defined above which 
can, be constructed with r nearest neighbour links in the lattice. 

The expansion, (34) can readily be generalized to the case of asymmetric 
interactions in different directions in the lattice. Thus for an 8.Q. lattice 
with a horizontal interaction J and a vertical interaction J’ and n 
horizontal and m vertical nearest neighbour links it is easy to show that 


Ziy! = 28 oosh Koos K'J"| 1+ s x p(r,s) tanh’ K tanhs K’ 
r=1s=1 


(35) 


where p(r,s) is the number of polygons with r horizontal and s vertical 
links. 
Fig. 7 


(a) (b) (c) (¢) 


Typical polygons in the zero field expansion. (a), (b), (c) contribute ; (d) does 
not contribute. 


‘Finally in terms of our previously defined functions and variables the 
relations (34), (35) can be written, 


Ay(ly2)=2" (+ ee [2+ eee |, v= a-ayaase. 
Il were) = Sol E Yolr, spor’ | : 


T, 8 


(36) 


the second relation applying to a quadratic lattice with equal numbers of 
horizontal and vertical links. 


3.4.2. The duality transformation in two-dimensional nets 

Using a unitary transformation on the matrices discussed in § 3.2.2 for 
the two-dimensional quadratic lattice, Kramers and Wannier (1941) 
established the relationship : 


2 Pelz) 
IL 2) acai t° Eee ee (ae) 
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They were thus able to establish a 1-1 correspondence between the partition 
function for a value z, and the partition function for a value v = (1 —z)/(1 a z). 
As T->0, z>0 and v1, and vice versa; knowing the partition function 
at low temperatures from z=0 to z=%e=(4/2—1)=(1—2e)/(1+2c) the 
values from z= to z=1 could automatically be deduced. Kramers and 


Fig. 8 


oP PP eeeees,. 


(b) 


Wannier thus concluded that if any singularity occurred for z < ze it would 
be repeated for z> Ze, and if there was only one singularity it must occur 
at z=%. Hence they conjectured that the Curie point of the quadratic 
lattice occurred at z=ze=(4/2—1). By evaluating the specific heat at 
this point for finite matrices of increasing magnitude, Kramers and 


Cooperative Phenomena in Crystals 179 


Wannier also conjectured that the specific heat was infinite at the Curie 
point. 

A topological interpretation, of this transformation, was given, by Onsager 
(see Wannier 1945) which gives a clearer understanding of its nature, and 
enables its application to be extended considerably. For any two-dimen- 
sional net we define a dual net by placing spins into every elementary 


Fig. 8 (continued) 


(d) 


Two dimensional nets and their duals (after Wannier). 


polygon, and drawing lines of interaction between them so that each old 
line is crossed by a new one. The duality relation is clearly a reciprocal 
one; the number of polygons in one net equals the number of vertices 
in the other, while the number of sides is the same. Examples of two- 
dimensional nets and their duals are given in fig. 8, and have been taken, 
from Wannier’s paper (1945). In cases (b) and (c) the dual net is 


180 C. Domb on 


topologically identical with the original; such nets may be referred to 
as ‘self dual’. 

For any Ising net it was pointed out in § 3.4.1 that the partition function 
can be expanded as a series of. terms, each corresponding to more over- 
turned spins and hence more deviations from the ground state. Thus for 
an, Ising net of N spins and / links in the absence of a magnetic field, 
we may write 

Zy=2exp(WVdI/kT) [: + > v(r)exp — (rd [kD) lt . (38) 
r=1 

The factor 2 arises because of the degeneracy of the lowest energy state 
and v(r) is the number of configurations corresponding to r broken links. 
It is not difficult to see that the number of configurations v(7) with r broken, 
links in a given, lattice is equal to the number of closed polygons p(r) of 
the type defined in § 3.4.1 for the dual lattice (see fig. 9). Thus if we use 
an asterisk to denote parameters of the dual lattice v(r)=p*(r). From 

(34) we have 

NM* 

Zn = 25" (cosh Ke 7 [ +> p(r) sane K* | oS ana 


r=1 


Fig. 9 


Equality of broken links on original lattice with polygons on dual lattice. 


Hence if we relate K* to K so that 
tanh K*=exp(—2K) ..... . (40) 
we obtain 
Zy*(K*) Zy(K) 
24* (cosh K*)“* a 2exp(WVK)’ 
Using the topological relations WM = *, N+N*= # +2, (41) can be 
put in, the more symmetrical form 


Zy(K) s Zy*(K*) 
24/2 (cosh 2K)42 — 2N*/2 (cosh 2K*)4?2 ° 
The relation, (42) means that if we know the partition function, of any 
net we can immediately write down, the corresponding partition function 


(41) 


(42) 


ee 
+ We shall drop the superscript J in Zy when this does not lead to confusion. 
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for the dual net. The relation (40) may be put more symmetrically in 
the form 

pinion tsinhok eles) eae, 4 (43) 
or in terms of z variables 


z*=v=(1—z)/(l+z); (l+z)(l+v)=2f. . . (44) 

The transformation, (37) is derived from (42) when, the dual net is identical 

with the original net (with the exception, of end effects) and this is the case 
for a quadratic net (fig. 8b). We find that 


ON/2 h 2K* N'|/2 


QN*/2 
and (37) follows readily (V*~N~ W/2) on taking the Nth root and 
substituting from (44). 

When a quadratic lattice is wrapped on a torus the relation v(r) = p*(r) 
is satisfied only for configurations which do not loop the torus (or loop 
it an even number of times) ; the duality relation is, therefore, not satisfied 
exactly. 

The duality relation, (42) can be established equally well for asymmetric 
nets using (35). Corresponding to (44) we will now have relations 


O= (be )/(l 2). v= (1=2)/(i-+2), . -. 2. . (46) 
and instead of (37) we find 
| TN Ca A z ENGL ye?) 60): Be ene (a7) 


(l+v)(1+2’) 
The relation (46) can no longer be solved for a given set of interaction 
constants to provide a unique relationship between 7'* and 7. However, 
if we interpret the transformation from z to v variables as converting order 
into disorder the Curie points correspond to the meeting points of the 
ordered and disordered regions in the (z,z’) plane; they will therefore be 
determined by 
g=(1—2'V(1+2'), 2’ =(1—z)/(1+2), _ | 
or (43) 
(l+z)(1+2’)=2 

Using methods similar to Kramers and Wannier, Potts (1952) established 
an inversion, transformation for a generalized interaction allowing more 
than two configurations for an individual spin. If the spins are allowed 
three configurations each 120° to one another, 


aS 


+The variable v will be used for the inversion transformation § 3.4.3. For 
the quadratic lattice the duality and inversion transformations are identical. 


P.M.S. O 
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and —J, is the interaction between parallel spins and —J, between 
spins oriented at 120° to each other, the Curie point is given by 


z=exp[—(Jy—J,)/kATcJ=(1+/3)*. . . - (49) 


This is a special case of a more general problem in which r kinds of units 
in equal abundance form a square lattice, and nearest neighbours interact 
with only two possible energies, — Jo for like units, and —J, for unlike 
units. For the general problem the Curie point is given by 


zo =exp[—(Jy—J3)/kTceJ=(1 +r). «© «. « (50) 


The particular case r=4 had been previously considered by Ashkin and 


Teller (1943). 
If the spins are allowed four configurations at right angles to one another 


| 


<—_-—> 


and —J is the interaction energy of parallel spins, +J of anti-parallel 
spins, and zero of perpendicular spins, the Curie temperature 7’, is given by 


Ze= exp (—23/kT'-)=(84+24/2)4, . .. «. (51) 


and is exactly half of the Curie temperature for the simple Ising model. 

These particular models do not correspond exactly to an Ising inter- 
action with a higher value of spin. However, since the duality transforma- 
tion, can be applied to them, it is rather likely that an exact solution can 
be found in closed form (as for the Ising model), and they might, therefore, 
provide a useful indication of the effect of a change in the type of inter- 
action, on the thermodynamic properties of the model. 


3.4.3. The star-triangle transformation 


The dual of a triangular lattice is a honeycomb lattice (fig. 8), and 
hence the relation (42) cannot be used directly to determine its Curie 
point. However, Onsager made use of the ‘ star-triangle’ transformation 
to relate the partition functions for the triangular and honeycomb 
lattices, and hence to determine the Curie points of both. The star- 
triangle transformation, replaces a star consisting of three spins interacting 
with a central spin, by a triangle of spins interacting with one another 
(see fig. 10); it can readily be established by using the relation (32). 

Let us consider the simplest case in which all the interactions are equal. 
Then summing over oy in, the partition function for the star, we obtain 

> exp K(o90, + o9¢_ + 0903) = 
a= +1 


> cosh? K(1 + ogo, tanh K)(1+ 090, tanh K)(1+ oo; tanh K) 
oo=+1 


=2cosh® K +2 sinh? K cosh K(o,0+ 0203 +0303). (52) 
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However, the partition function for a triangle of spins with interaction 
K* is given, by 
exp K+(o409 + 0203 + 0304) = 
cosh? K*(1 + 0,0, tanh K+)(1 + 0,0, tanh K+)(1 +030, tanh K+) 
= (cosh? K* + sinh? K+) + (cosh? K+ sinh K++ cosh K+ sinh? K+) 


X (G40g+0903+030;).. . . (53) 


Fig. 10 


The star-triangle transformation. 


Comparing (52) and (53), we see that if we chose K+ to satisfy 
cosh K+ sinh K+ (cosh K++ sinh K+) 
(cosh? K++ sinh’ K+) 


the interactions of the star and triangle will be related to each other by 
the simple factor 


=tanh? K . . (54) 


f=2 cosh’ K/(cosh® K++ sinh? K+). 
If we now apply the star-triangle transformation, to the lattice points 
represented by open circles in, (fig. 11), we easily deduce that 
Zoy"(K) a ein te) . . . . . . (55) 


where the superscripts H and 7’ refer to the honeycomb and triangular 


lattices. 
Now from (42) we have the duality relationship 


24/2 (cosh 2K*)34/2 2" (cosh 2K) 4/2 ° 7 ee 
Combining (55) and (56) we obtain a relation between Zy*(K*) and 
Z,"(K+), to which the argument of the previous section can be applied. 
The relation (54) between the variables can be more simply expressed as 
exp 4K+=2cosh 2K —1, \ Pe Br) 
at? = 3/(l—-z+2*), 
02 
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and hence from (44) the reciprocal relation between 2+ and 2* is 


semicon, 8 
[exp (4K+) —1][exp (4K *) -1]=4. 

Following Wannier (1950) we shall refer to the transformation between 
variables such as z* and z+ as an ‘inversion’ transformation. For the 
quadratic lattice the duality transformation is thus identical with the 
inversion transformation. Using (43) and (45) it can be put in the 
alternative form 

Zyl KY oe DE 
(2sinh2K)*2 (2sinh2K*)42" 
The Curie point of the triangular net is given by 2*=z*, or 
(1+2¢?)(1—32e2)=0, ze=1/4/3. 


(59) 


Fig. 11 


The honeycomb lattice. Application of the star-triangle transformation to 
the spins represented by open circles and their three neighbours leads 
to the triangular lattice (broken lines). 


The inversion transformation for the triangular lattice follows from (55) 
and (56); Wannier noted that this transformation was of exactly the 
same form as (59) in the variables K+ and K*. The form corresponding 
to (37) is 


A(l,2)= A(1,v), v?=(1—2?)/(1+322) . . . (60) 


1+3v? 
and this form is useful in the derivation of series expansions. 

For the honeycomb lattice the corresponding relationship from (55) 
and (56) is between Z,y"(K) and Z,,y"(K+*) where 2+* =(1—2+)/(1+2+). 
The Curie point is given by z=zt* so that (1+2¢%)(1—4z2e+2¢2) =0, 
%e=2— 4/3. The relation between z and zt* is 

gtk — V/(1—2+2%)—+/z 
V/(L—2+2*) + 4/2’ So ag OES 
(cosh 2K+* —1)(cosh 2K —1)=1. 
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Again the relation (59) is valid in the variables K and K+*, and the trans- 
formation can be written 


A(1,2z)=2[1 +p it (=a) |” A(1, 2), 


1—v+v? 
Ree a) Ve 72) 
V(L-2+2)+ Vz" 


It will be noted that the equations determining the Curie points of the 
quadratic and honeycomb lattices have identical solutions for z and 2-1, 
whereas the triangular lattice has only the one solution ze=1/1/3. This 
corresponds to the fact that perfectly ordered antiferromagnetic arrange- 
ments exist for the first two lattices but not for the third (§ 3.5.4). 

Duality transformations can, be derived in a similar manner when, the 
interactions are different for different directions in the lattice. The 
detailed calculations are somewhat complicated by the methods of this 
section, and follow more readily from the general approach of Fisher to 
be described in (§ 3.4.5). We shall quote the results for reference (Houtappel 
1950, Potts 1951). 

For a triangular lattice with ferromagnetic interactions J, J’, J” the 
Curie point is given by 


pea eye le eee Gt), x (63) 


When one of the interactions becomes zero we recover the formula (44) 
for a rectangular lattice. More generally when interactions can, be ferro- 
magnetic or antiferromagnetic the equation determining the Curie point is 


sinh 2K sinh 2K’ +sinh 2K’ sinh 2K" +sinh2K"’sinh2K=1. (64) 


This equation is valid when the largest two of J, J’, J” are positive, i.e. 
J+J’>0, J’+ J” >0, J”+J>0. Other cases can be dealt with since 
any two of J, J’, J” can change sign without affecting the Curie point, 
and (64) can then be used for the new values if they satisfy the condition 
of validity. There is no Curie point if and only if one or three of the J’s 
are negative and the two weakest J’s are exactly equal; in these cases 
no ordered lowest energy state occurs. 

For the honeycomb lattice with different interactions in different 
directions the corresponding equation is 


sinh 2|K|+sinh 2|K’|+sinh 2|K”|= 
sinb 2|K|sinn 2|K’|sinh2|K"|. =. . (65) 
This lattice always has a Curie point for any values of J, J’, J". 


3.4.4. The decoration transformation 

A transformation similar to the star-triangle transformation was 
suggested by Syozi (1951) and developed by Naya (1954). This applies 
to a single central spin connected to two spins, and replaces them (including 
the interaction with an external magnetic field) by an interaction along 
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a single bond between, the two atoms, and a change in magnetic field inter- 
action of these atoms (see fig. 12a); the central spin will be referred to as 
a decorating spin. 
Taking the simple case shown in fig. 12a, and summing over op in the 
partition function for the three spins we obtain 
> exp (Koyo, + Koyog+ Loy), (L=mH/kT) 
o=+1 
= > (cosh K +o 90, sinh K)(cosh K + ogo, sinh K) 
%=+1 
x (cosh L + op sinh L) 
= 2cosh? K cosh D + 20,0, sinh? K cosh L + 2(a, + 09) 
x cosh K sinh K sinh L. Ste Se GG) 


Fig. 12 


(a) 


General 
Arbitrary 
Interactions 


The decoration transformation. (a) Simple transformation for atoms of spin 
1/2; (6) Central spin replaced by an arbitrary interacting system. 


But the partition function for a pair of spins with interaction K and 
magnetic interaction L is given by 
exp (Ko,0, + Lo, + Loy) = (cosh K + 0,0, sinh K)(cosh L +, sinh L) 
x (cosh [+ o, sinh L) 
= cosh K cosh? + o,0,[cosh? K sinh L + sinh? K cosh L] 

+(o,+0,)(coshK+sinhK)coshZsinhZ. . . . . . . (67) 
Hence K and L can be chosen so that the two expressions (66) and (67) 
are equivalent except for a simple factor. The appropriate values of 
K and L satisfy. 


tanh? K = tanh K + tanh? L, 6 
tanh K tanh L=tanh Z + tanh K tanh L. } ee, 
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Using this transformation Naya derived a relation between the partition 
functions for the magnetic honeycomb lattice with parameters K, L, 
and the decorated honeycomb lattice with parameters K, L as shown in 
fig. 13. (In the decorated lattice spins represented by open circles are 
non-magnetic and do not interact with an external field.) The star- 
triangle transformation was now used to remove the non-magnetic spins 
and the kagomé lattice (fig. 13 6) was derived. He thus obtained a relation 
between the partition functions of the honeycomb and kagomé lattices 
in the presence of a magnetic field. The relation may be symbolically 
characterized as follows: 


Pe) Ze OK, LZ AtKe, Ly). . . « (69) 


The relation (69) can be used to derive the Curie point of the kagomé 
lattice. Taking L =0 we find from (68) that L =0, and the relation between 
K and K reduces to - 

HTN CA 2a EW el ech Le ls =) Vee weno ves ea OC K8)) 


The decorated honeycomb lattice and the kagomé lattice. Application of the 
star-triangle transformation to the spins represented by open circles 
(a) yields the kagomé lattice, (0). 


Using the relation (49) between z and 2+, we find for z+ and z 

Derg Oras se ES ge (TY) 
and since the Curie point of the honeycomb lattice is given by 7 =2—+/3, 
the Curie point of the kagomé lattice is given by 


9 1 1/2 
(e5) 
(Syozi 1951). 


The kagomé lattice has coordination number 4, and all points and. 
bonds are equivalent, but its Curie temperature is lower than that of the 
quadratic lattice. 
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3.4.5. Generalized transformations 

Recently Fisher (1959 a) has put forward an alternative approach to the 
decoration and star triangle transformations which appreciably simplifies 
the algebraic manipulations, and enables the transformations to be 
generalized considerably. Let us first consider the decoration transforma- 
tion but replacing the decorating spin by any arbitrary system interacting 
with the lattice spins (fig. 126). The energy levels of this system can then 
be written in the form H,(c,,0,; H), and the partition function (summing 
over the decorating system) by 


Hoy, 02; H)= DYexp[—Ho,,0,; Hk]. —. «. (72) 


This is a function of the spin variables o,, c, taking on four values corres- 
ponding to the four possible combinations (+, +)(+,—)(—,+)(—., —). 
We wish to replace this function by an interaction along the original bond 
K, together with a change in the magnetic interactions of the lattice 
spins L,, L,, ie. by 

fexp (Ko,o,+-L,0,4+ 1,0.) - 9.) so ee 


where f is a factor independent of the spin variables. There are four free 
variables in (73) and four independent conditions to be satisfied in equating 
to (72). In fact it is easy to see that they will be satisfied by taking 


(74) 


i t= $i p_ pb, p_,, 
where #,,=%(+1, +1; H) and so on. The example discussed in the 
previous section is a special case, and K and L are explicitly determined 
by (74), 
cosh (2K + L) cosh (2K —L) 
cosh? L ; 

exp (4Z) = cosh (2K +L) cosh (2K —L), 

f* = cosh (2K +L) cosh (2K —L) cosh? L. 


exp (4K) = 
(75) 


But corresponding formulae could readily be written down for different 
decorating systems, e.g. a spin of magnitude s, or a group of spins. 

The partition function for a decorated lattice is reduced to a summation 
over the original undecorated lattice 


DK, Ly=f* ZK, BE), se te) 
where 4) is the number of decorated bonds. It may be noted that the 
factor f is effectively the partition function of the decorating system 
averaged over all possible states of the basic lattice spins to which it is 


coupled. Consequently it will only give rise to critical behaviour if the 
individual decorating systems exhibit critical behaviour; otherwise the 
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singularities of the decorated lattice will be simply related to the singulari- 
ties of the original lattice. 


Using the same approach towards the star-triangle transformation, 
and replacing the central spin by an arbitrary interacting system, we must 
consider the energy levels H,(o1, 0,03; H). In general there are not enough 
free parameters in the partition function for the triangle to enable us 
to establish an identity valid for all o,, o5, og. However, if we restrict the 
central interacting system so that the partition function 


P(o1, F2, 3) = 29) eB iket og os) RT| 4. a 3(77) 
is Invariant under the operation of total spin inversion 
Chia Ohp 09> — 09, 03-7 — Og, 
we can establish such an identity. Writing 
W+,+,+)=$—, —-, —)=po, 
Ws — +)=W-, +, —)=po he 
Ure teat) UN pas hh) =, 
we must equate (77) with the partition function for a triangle of spins 
FOxpi Ke 0,06, 050, + Kaie,a,).- . . . «= (79) 
This can, be satisfied identically if we take 


exp (4K,+) =d5 ie exp (4K 4+) =p Yo 


o's Pe’ en 80) 
exp (4K37) = Se , FA = poprporbs: 


Again, the partition function, of a ‘starred’ lattice is related to the partition 
function, of the basic ‘triangle’ lattice by 

ZS(Ky, Ky, Ks)=f AK y*, Ky*, Ks"), 
where NV, is the number of star vertices. 

The restriction of invariance in the central interacting system means 
that there is no preferred direction for this system, and effectively restricts 
the spins of this system to having no magnetic moment. As an example 
of the application of (80) we can write down the general form of the star- 
triangle transformation, for a single central spin with different interactions 
in different directions discussed at the end of § 3.4.3: 
cosh (K +.K’+ K")cosh (—K+K'+K") 


De ee ee 
aa cosh (K —K’+ K") cosh (K + K’— Kk") 
4R'+) = cosh (K + K’+ K") cosh (K — K'+ K") 
eEP a cosh (— K +.K’+ K") cosh (K + K’—K") ’ (81) 
/ " h K Ke K 
expiikia)= cosh (K + K’+ 4") cosh (kK + ) 


cosh (— K + K’+ K") cosh (K —K'+ Kk") 
ft=2*cosh (K+ K'+K")cosh(—K+K'+K", 
cosh (K —K’+ K")cosh (K+ K’— K"). 
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Similarly for a central spin of magnitude s interacting with three standard 
spins of magnitude } the transformation energy parameter 1s given by 


exp (4K+)={1+2cosh [(2s + 1)K/s]}/[1+2cosh(K/s)] . . (82) 


where the interaction energy between, the central and outer spins is 
normalized so that its maximum and minimum values are + J(K=J/k7). 
By this means we can relate the partition function of the alternate 
honeycomb lattice (fig. 14) to that of the triangular lattice. 


Fig. 14 


Honeycomb lattice half the spins having s=1/2, and the other half having 
arbitrary s. 


3.5. Exact Solutions for Two-dimensional Ising Models 


The most important landmark in the field of cooperative phenomena 
was the exact solution in 1944 by Onsager of the two-dimensional Ising 
model in, the absence of a magnetic field. Onsager’s method involved the 
theory of group representations, and a substantial knowledge of algebra 
is required to follow the detailed work. A few years later Kaufman (1949) 
substantially simplified the development by the use of spinor analysis, and 
Kaufman, and Onsager (1949) showed how the degree of order could be 
calculated between any two lattice sites. In 1952 Kac and Ward intro- 
duced a more elementary approach to the problem, which can be understood 
without detailed algebraic knowledge; they were able to reproduce all 
the results which had been derived previously. It is not our intention to 
provide the mathematical details of the various derivations, since they 
are Clearly set out in the original papers, and have been further discussed 
ina review paper (Newell and Montroll 1953). However, we shall endeavour 
to summarize briefly the various methods of approach, and to reproduce 


fully the conclusions regarding the mathematical and physical properties 
of the resulting solution. 
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3.5.1. Onsager’s solution 


Onsager considered a two dimensional lattice each layer of which 
contained spins connected to form a circle (fig. 15). The theory of §3.2 
applies, and since successive layers are connected to form a ring, the 
complete lattice is wrapped on a torus. The interaction U(c,o+) between 
successive layers, which gives rise to a matrix V of order 2”, can, be divided 
into two parts, a portion between, successive layers, and a portion within 
the added layer. It can be written symbolically as 


U(c,o7)=—J > o0;*-J' > o6,,, (C4n=o;), - ~ (88) 
y i=1 


t=1 
where the o; are variables which can take the values +1, as in §3.3. 


Fig. 15 


~~ = — = 
- — 
oo CF Ses 


Each added layer consists of a circle of n spins. The added interaction can 
be divided into two parts (i) a portion between the final two layers 
@ and O, (ii) a portion within the final layer O. 


Correspondingly V can be written as a product 


n 


V= J[ exp K(o,0;*) [J exp K'(o,0;,,). eave + (84) 
i pa 


i=1 
The exp K(c,o;+) can be expressed in terms of operators with simple 
properties. When only one particular 7 is considered exp K(g,o;+) can be 
represented by a 2 x 2 matrix. 


K —K 
(C= *) =eXI +e-KC, | e=(; ah ee (85) 


The operator C; changes o; into —o; and satisfies C?=1; thus using (43), 
(85) may be written in the alternative form (2sinh2K WZ exp (K*C,). 
When, n independent o; enter, each C; can be represented by an operator 
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in 2” dimensional space, and C;, C; commute. If we represent any eigen- 
vector by |o,,..., oj...) we can write 


Ciloy, eas Ot a oe lois ty Oy Geen Os vd So) 
We can define similarly an operator s,; which does not change the sign of 
o;, so that it is represented by 
(0 -1) 
0-1)" 


Slo) «<2 Ops = 21 Oy ) ONES Phas Crees Gare . . . (87) 


Corresponding to (86), we have 


In terms of the operators C;, s;, V can be written in the form 
(2 sinh 2K)"?V,V., 
where 
n 
Vi= II exp (K*C;), 


i=1 


n 
V.= I] exp (K’s;s;,,). 


i=1 


(88) 


Alternatively if we wish to retain the symmetric property we can take V 
in the form (2 sinh 2K)"/?V,!?V,V,"?. The commutation properties of the 
C,, s; are as follows: 


s?=CZ=1, sC;= —Cs;, 
The set of operators (1, s;, C;, s,C;) form a simple quaternion algebra, and 
hence the algebra of n such sets is a direct product of n simple quaternion 
algebras. 
Defining operators A, and A, by 


ee Sire 


i=1 


n 
A,= > $jSii4, 
i 


v 


i? 


(90) 


Onsager proceeded to form the commutator [A,A,], and then the 
commutators of this with A, and A,; he continued to form commutators 
of commutators until he obtained a complete set of operators such that 
the commutator of any two operators is a linear combination of members 
of the set. Such a set of operators is called a Lie algebra. The Lie algebra 
generated was of order only (3n—1) and possessed certain characteristic 
symmetries. Making use of these symmetries Onsager was able to reduce 
V,V, to a product of commuting matrices which could te diagonalized 
separately (for details see Onsager 1944, Newell and Montroll 1953). 

If the same procedure is attempted for a two-dimensional lattice in 
the presence of a magnetic field, or for a three-dimensional lattice, the 
Lie algebra generated is much larger, and the simplifying features are 


Cooperative Phenomena in Crystals 193 


no longer present. No progress has yet been made in the algebraic 
approach to these more general problems. 


The basic properties of Onsager’s solution may be summarized as 
follows : 


(a) Mathematical properties of the solution 
(i) The largest ergenvalue Amax is given, by the formula, 
In Amax — $7 In (2 sinh 2K) =4(y, +y3+...+Yen-1)- . (91) 
The y, are given by 
cosh y, = cosh 2K * cosh 2K’ — sinh 2K * sinh 2K’ cos rz/n 
= coth 2K cosh 2K’ — cosech 2K sinh 2K’ cosra/n (92) 
=coth 2K cosh2K—cosrz/n, (K=K’). 


It will be seen that the y, are invariant under the duality transformation 
K’+K*, and hence so is the sum in (91); they are shown graphically in 
tee eG (ACK), 


Fig. 16 


ie) A me Re) 4 5 6 ot 8 
K 


Dependence of y, on temperature (after Kaufman). For 7>T7 (kK <K,), 
Yiz—Y%- For T<T.(K> Ke), yi=Yo- 


In the limit when, the size of the crystal becomes large the series can 
be replaced by an integral, and the partition function per spin is given by 
; vee 
Woe tin ainh K)= =| Par (98) 
0 


This can be transformed into a double integral which is symmetric in 
K and K' [using the identity 


2 . 
| In (2 cosh a — 2 cos w) dw = 272], 
0 
In (A/2) = a["{" In (cosh 2K cosh 2K’ — sinh 2K cos w 
oJ 0 
—sinh 2K’ cosw’) dw dw’. een.c eo ae + on (D4) 
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The various thermodynamic functions can be derived in the standard 
manner from A, and many of them are expressible in terms of elliptic 
integrals; the results are summarized in Appendix II. 

(ii) The distribution of eigenvalues. Onsager also obtained the remaining 
eigenvalues of V, and these enable us to understand the ordering properties 
of the assembly (§3.3.1). The eigenvalues fall into two sets (Kaufman 
1949) 

IndA,t+— 4nIn (2sinh 2K)=4(+y, ty3---+Yon—1)s \ _ (95) 
Ind; — 4nIn (2sinh 2K)=4(+yotya---+ Yen): 
Here any combination giving an even number of minus signs corresponds 
to an eigenvalue; there are thus 2” eigenvalues as expected. 

It will be seen from fig. 16 that the y, increase steadily from y, to y, 
(=K'+K*), and for r>n they satisfy the relation y,=y2,_,; for large n 
they form a closely spaced spectrum between 2|K’—K*| and 2(K’+K*). 
Yo is exceptional, however, since it changes sign at the critical point K’ = K*. 


Fig. 17 


Schematic representation of the distribution of eigenvalues. Shaded portion 
corresponds to a continuum as > 00. 


For large ”, Yor41~ Yor (r#0) and for K’> Ko, ie. T<To, yy~% yo; thus 
2a 
the two sums in (95) are different approximations to} | y(w) dw, and the 


largest eigenvalues A,+ and A,~ are very nearly equal. (The difference is 
exponentially small, all correction terms in the Euler—Maclaurin formula 
vanishing; see Goodwin 1949.) The next eigenvalue differs from this 
value by a factor of approximately exp4(K’—K*). However, when 
K’< Ko, i.e. T>T'c, yy~% —yo. Thus the largest eigenvalue A,+ is greater 
than A,~ by a factor of approximately exp [2(K*—K’)]. 

We may summarize the distribution of the eigenvalues as follows: 
T<Tc, A,+=A,~ > continuum of limit \,+ exp [—4(K’ — K*)], 
T=T., A,+=A,-=limit of continuum, 

T>To, A+ >A,~=limit of continuum (\,-=A,+ exp [— 2(K* — K’)]). 
This is illustrated diagramatically in fig. 17. 
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Following the discussion of order in § 3.3 we see that for 7’ < 7. there is 
long-range order, arising from the degeneracy of the largest eigenvalue 
(but this occurs only in the limit as n+0o). For 7'>7'c there is no long- 
range order, and the mean distance to which order persists is given roughly 
from (26) by the inverse of In (A,+/A,~) or 1/[2(K* — K’)]. 

(iil) Change in sign of the interaction. If we replace J’ by —J’ the 
interaction, becomes antiferromagnetic in one direction; the solution, (95) 
remains valid but a detailed analysis shows a difference between odd and 
even values of n. We may now write 


cosh y,= cosh 2K * cosh 2K’ — sinh 2K * sinh 2K’ cos(m—rz/n), 
(96) 


so that each y, is replaced by y,_,.. When is even),+ and A,~ are unchanged 
as we should expect. When n is odd, however, there is a misfit seam; 
odd y’s are transformed into even y’s, except for y,(=K’+K*) which is 
transformed into —y(=—K’'+*). Thus ),* is transformed into A,~ 
with the exception of the change of sign of y,». We may, therefore, use 
the discussion, of (ii) to conclude that when 7' > 7. the partition function is 
unchanged, but when 7'< 7, there is an additional factor 


exp [—2(K’— K*)]; 


this will be interpreted as a ‘boundary tension’. However, the largest 
eigenvalue becomes multiply degenerate for 7’<7'c, as we should expect 
from the degeneracy in position of the misfit seam. 


(b) Physical properties of the solution 


(i) The specific heat manifests a striking singularity at the Curie point. 
It is logarithmically infinite on both sides of the Curie temperature, and 
is drawn, in fig. 18 for the symmetric net. For temperatures near to 7’. 
the following approximation, may be used: 


Cy |/nmk ~ (2/)(kn cot 7/8)2(K, — 1— 7/4), 
K,=In (21?/|K — Ke). 


It should be noted that the specific heat curve differs substantially from 
any curves which have been observed experimentally, particularly in 
regard to the magnitude of the ‘tail’ above the Curie temperature. This 
can, be seen, more clearly from the critical value of the entropy, which 
shows that the entropy change above 7’, (=0-387) is greater than, that 
below 7’, (= 0-306). Itis important to assess whether this is a characteristic 
of the Ising interaction, or merely dependent upon the two-dimensional 
character of the model; we shall take this point up again in §4.7 in 
discussing the properties of the three-dimensional model. We may 
incidentally draw attention to a recent remarkable experiment by 
Fairbank et al. (1958) on the anomaly in liquid He which has shown clearly 
that it is logarithmic on both sides of the critical temperature. 


(97) 
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For an asymmetric net the logarithmic nature of the specific heat 
anomaly is unchanged. As one of the interactions J’ becomes small the 
Curie temperature moves to zero, but this movement is very slow; if the 
interaction J’ is J/v then the reduction, in 7’; is of order Inv for large v. 
The specific heat curve for an asymmetric net in which J’=J/100 was 
calculated by Onsager, and is shown, in fig. 18 for comparison with the 
symmetric net, and with the specific heat of a one dimensional chain 
(J’=0). The abcissa is 2k7'/(J + J’) so that the areas under all the curves 
are the same. It will be seen that the large asymmetry reduces the Curie 
point by a factor of 2 only. Also there is an essential difference below the 
Curie temperature between a two-dimensional net with one small inter- 
action, and a one dimensional chain; as pointed out in § 3.2.1 the latter 
can be regarded as an assembly whose Curie point has moved to 7’'=0. 


Fig. 18 


Specific heat of a rectangular lattice. -—----—--+-- J’=J (symmetric) 
——_J'=J/100; ------- ‘=0 (linear chain) (after Onsager). 


(ii) Other thermodynamic quantities. The internal energy EF has a 
vertical tangent on, both sides of the Curie point and is shown in fig. 19 
(for the symmetric net). Near the Curie point it is of the form 

A,+B,(T—T-)In|T —T-I. 
The free energy F involves an integral of the internal energy and therefore 
has infinite curvature at the Curie point. The entropy can readily be 
obtained from the relation S=(H#/7)—F. Its behaviour at the Curie 
point is similar to that of theenergy. Critical valuesare givenin Appendix II. 

(ili) The boundary tension. We have seen that when, the sign of the 
interaction, energy J’ is changed, if n is odd there is a misfit seam which 


divides regions of opposite order. The partition function (per spin) of 
such a boundary is 


A, =exp[—2(K’—K*)]=2'(1+2)/(l—z). . . . (98) 
The corresponding boundary free energy is 
o =2J'—kT Incoth(J/kT); 
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it vanishes at the critical point and is shown in fig. 19 for a symmetric 
net (J=J"). It is thus significant that the critical points are identical for 
the boundary and for the bulk substance. Corresponding boundary 
functions A,, o in the perpendicular direction are obtained by interchanging 
J and J’. 

It was pointed out by Onsager and Temperley (Temperley 1952) that 
the formula (98) is exactly the partition function which would be obtained 
if only a particular class of boundaries (forward boundaries in the J’ 
direction) were considered. Temperley (1952) attempted to use this result 
to draw some conclusions about polymer configurations on a quadratic 
lattice, but Fisher (1958) showed that these considerations are invalid. 


Fig. 19 


2:0J 


1-OJ 


Energy and boundary tension of a quadratic lattice. 
Hf has imiinite slope ati 73 +--+ «2-** > +: o is zero for T> Ty. 


(iv) The specific heat of a finite crystal. The specific heat of a finite 
erystal of n rows can, be computed from eqn. (91). It is finite at all tempera- 
tures, since the partition function is an analytic function of the temperature, 
but a maximum occurs near the transition point of the infinite crystal ; 
this maximum increases and sharpens with increasing size of the crystal. 
Onsager showed that for large m the maximum specific heat is given 
approximately by 

Cy/nmk =~ 0-4945 Inn + 0:1879, SPP Te pea 199) 
and it occurs at a temperature differing by 0(n~*Inn) from 7’c for the 
infinite lattice. The formula extrapolated by Kramers and Wannier 
(1941) from calculations for finite crystals with n=1, 2, 3, 4, 5, 6 is quite 
close to (99). 


P.M.S. 
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(v) The mean ordered length of a finite crystal. For an, infinite crystal 
below the Curie temperature there is long-range order, and the majority 
of spins will point either in direction 1 or in direction 2. However, for a 
strip of finite width this is no longer the case, and for a long strip equal 
numbers of spins will point in each direction on the average. We should 
nevertheless expect them to fall into blocks of opposite order, and the 
formulae for the boundary tension enable us to estimate the mean length 
of such an ordered region. The partition function of a boundary of finite 
width n between, regions of opposite order is approximately 4,” and this 
represents, by Boltzmann’s principle, the probability of such a boundary 
occurring at a given, section of the strip at temperature 7’. Hence the 
mean, length between such boundary occurrences is given from the theory 
of recurrent events in probability (Feller 1950) by 


lj=A,—” = exp [2n(K —‘K’*)]=[exp(2K)tanh K’)". . (100) 
Clearly when, n is large /,>n, and the estimate (100) is then significant. 


3.5.2. The development of Kaufman and Onsager 


In 1949 Bruria Kaufman showed that the algebra required for the 
determination of the eigenvalues of V was that of the spin representations 
of the orthogonal group, the general theory of which had been formulated 
previously by Brauer and Weyl (1935). The matrix V is a 2”-dimensional 
spin representative of a 2n-dimensional rotation, (i.e. orthogonal matrix), 
the angles of rotation being complex. The problem is thus reduced to the 
diagonalization of a 2n-dimensional orthogonal matrix, and considerations 
of symmetry resolve this into a product of n plane rotations. The eigen- 
values are again, found from the solution of » quadratic equations. In 
addition the matrix of eigenvectors which diagonalizes V is determined, 
and this can be used to calculate the degree of order at various distances. 

Although the procedure may be standard, it is difficult for anyone without 
specialized algebraic knowledge to follow. We shall, therefore, again 
summarize the main conclusions. 


(i) The exact partition function for a torus. The complete set of eigen- 
values has already been given in (95). Hence the exact partition function 
for an m xn lattice wrapped on a torus can, be written down, from (2) 


n n 
Zm,n= 9(2 sinh 2K yr [] 2cosh (my.,_,/2)+ [] 2sinh (myp,_4/2) 
1 r=1 


T= 


+ [] 2cosh(my,,/2)+ T] 2sinh (myr/2) | : 2 ee 
r=1 


r=1 


For large n we have seen in §3.5.1 that yo..1> Yo, (r 0), 
Yory1 for T<T., and yy~—y, for T>T-. 


Therefore the two products containing cosh (my/2) as factors are approxi- 
mately equal (except for temperatures very near 7'c) while the two products 


a 
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with sinh (my/2) are approximately equal for 7<T'. but cancel out 
approximately for T’> 7. Thus we have 


Zn,n@ (2sinh 2K)? TT 2 cosh (myo,_,/2) 
r=1 


x [2 + |] tanh (my -s12) | (Peis); - (102) 
r=1 
Zn, n@ (28inh 2K yn" TT Qoosh (myo,4/2), (> Te). 
r=1 


Also for large m, 2 cosh (my>,_/2), exp (My2,_,/2), and the product of the 
tanh (my>»,_;/2) is equal to 1 (except in the immediate neighbourhood of 7). 


Hence 
OD \xce etl, URI Bh 
Los sat oe eee fee TOS 
ees { igs EL 6. } ( ) 


We thus see that the approximation, of considering only the largest eigen- 
value is a good one as long as m and n are large; the factor 2 arising for 
7 <T corresponds to the twofold degeneracy in Amax in the largest 
eigenvalue approximation. 

(ii) The degree of order. It was shown in § 3.3.1 that a knowledge of the 
complete set of eigenvalues and eigenvectors of V enabled the degree of 
order to be determined between any two lattice points. This can be 
conveniently expressed in terms of the matrices s, defined in (87). Each 
s, is a diagonal matrix of order 2” with +1 for all of the 2”—-1 components 
in which o,=1, and with —1 for all of the components in which o,= —1. 
Thus the trace of s, V” is obtained by multiplying all diagonal components 
of V™ with o,= +1 by +1, and all diagonal components with o,=—1 
by —1, and therefore 


Trace (s,V™) 
(is eee, 


Ae 
dat 


Similarly the correlation between, spin v in layer 1 and spin p in layer h is 
obtained from (23) by summing the product of the a8 component of V” 
with the Bu component of W"~" over all configurations « with a given 
value of spin v and f with a given, value of spin p. It is easy to see that 
trace (s,W’s, V"-") gives the factors +1 according to the appropriate 
values of the spins o,, and o,,, and hence that 


Trace (s,V’s,W"—") 


Qn 
Je 
v1 
Kaufman and Onsager applied their general theory of the diagonalization 
of V to the evaluation of these averages. The details of the algebraic 
manipulations are quite involved, and the resulting correlations are 
expressed as determinants of functions which can be reduced to elliptic 
integrals in the symmetric case K=K’. The nearest neighbour correlation 


(104 


{04,0 np) = (105) 


P2 
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(c,0,) is equivalent to the energy ; and as the distance apart becomes very 
large the correlation (0,0, ) approaches the long-range order 7 ,,, which has 
been evaluated exactly, and will be discussed shortly. All the intermediate 
correlations lie between these two bounding curves. Kaufman and 
Onsager evaluated exactly the correlations (0103) (0,022) (0193 ) {011923 ) 
and these are shown, graphically in fig. 20. At the Curie point it seems 
that all the correlation functions have a singularity similar to the energy 
with an, infinite slope. 

For the correlations (o,¢;,) along a given row a general formula can 
be obtained for all & in the form 


(—1)*(o,0;,4,)=cosh? K*A,+sinh?K*A_,, . . (106) 
+ k k 


Fig. 20 


kT/2J Tel 


Correlations for small distances as functions of temperature. (qa) (0102), 
(6) <o41092>, (Cc) 0403>, (d) <oy,023>, (€) <oyo,>. All have infinite 
slope at 7’, (after Kaufman and Onsager). 


Where A, and A_; are Toeplitz determinants of order k whose elements 
can, be reduced to elliptic integrals. At the Curie point the following 
approximations are reasonable : 


2 4-1 T(s+1)I(s4+1) 
A,|=- al a id A'S a 
ele UT Ter prer eh) 

2*=1 1(s+1)0(s+1) 
A_,|= = pat at ed Kt Gal 2 
Mole 5 OU perprer® 
Hence for large & the correlation at the critical point ~ A/k”* and tends 
to zero very slowly as k>oo. 


(107) 


tIn the paper by Kaufman and Onsager the products are taken to s=k b 
Dr. M. E. Fisher has pointed out that this is incorrect. tS i 
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Onsager (private communication) has also derived an expression for 
the order along a main diagonal for the general anisotropic case which is 
much simpler than (106). In fact if a, are defined by 


| sinh 2K sinh 2K’ — exp (iw) 


sinh 2K sinh 2K’ —exp(—tw) |= , =, #OmP ire), _ - (108) 


T=— © 


so that. the a, are Fourier coefficients of this function, the order along a 
main diagonal is given by a single Toeplitz determinant 


Gog 1, © Og Aya 
Os Cain) A. Ay_2 

D= A_» a_y Qo elie ‘eo 16 Ap_g . . . . (109) 
Q_(4—1) A_(7,_2) Yo fo Ay) 


For this determinant the first equation of (107) is exact so that the correla- 
tions at the critical point along a diagonal ~ 0-645/k"4 and the asymptotic 
form is approached very rapidly. 

In the previous section we derived an estimate of the mean range over 
which order persists above the Curie point in the form 1/[2(K* — K’)]; this 
is valid in the direction, of interaction J, and the corresponding expression 
in direction J’ is 1/[2(K*—K’)]. In intermediate directions of angle ¢ 
Onsager in, his first paper conjectured the implicit formula, 


cosh 2K cosh 2K’ — sinh 2K cosh (8 sin d) — sinh 2K’ cosh (8 cos ) = 0, 
. (110) 


for the mean range 1/f of the order. This would imply that, with K=K’, 
for temperatures near the critical point the order spreads out isotropically, 
whereas for high temperatures it spreads out along the links of the lattice 
(fig. 21). The exact computations available so far confirm this general 
behaviour. Knowledge of the mean range of order as a function of 
temperature provides information on the temperature variation of the 
correlation, functions, and this has been, used by Fisher (1959 b) to derive 
an, expression for the high temperature susceptibility near the Curie 
point. : 

(iii) The long-range order. ‘There are two independent approaches to 
the question of long-range order. The first considers the limit of the 
correlation function (c,0;,) defined above as k>oo. The second evaluates 
the spontaneous magnetization, by a small perturbation of the no-field 
solution, for 7’< 7c; the spontaneous magnetization, is equal to the long- 
range parameter F (defined in § 3.3) as will be shown in § 3.6.4. 

The exact formula for the long-range order was first given by Onsager 
(Florence Conference on Statistical Mechanics 1949) although he provided 
no details of his calculation. The result can be derived from (109) or 
corresponding determinants obtained by Potts and Ward (see §3.5.3) 
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using the theory of Toeplitz forms}. The only published derivation, (using 
the perturbation method) is due to Yang (1952) (see also Chang 1952). 
The result can be put in the remarkably simple form (Potts 1952): 


R§ =1-sinh? 2K’*/sinh? 2K \ Lt (o,,)=R. (LL) 
= 1— 162?2'2/(1 —22)2(1—z'2)? J pro 


R differs only slightly from unity except in the immediate neighbourhood 
of the Curie point; for the symmetric lattice its value near the 
Curie point is ~ 1-22241(1—7'/7.)"*. The spontaneous magnetization is 
represented graphically in fig. 32. 


Fig. 21 


(a) (b) 


Anisotropy in the propagation of order. At high temperatures order spreads 
out along the links of the lattice with little correlation and the contours 
of constant order are as shown in (a). Near the Curie point the correla- 
tion becomes very marked, and the contours are circles, as in (6). 


3.5.3. The combinatorial method of Kac and Ward 


The method of Kac and Ward is elementary in the sense that no detailed 
algebraic knowledge is required other than properties of determinants. 
It undoubtedly provides the simplest derivation of the partition function 
and other properties of two-dimensional Ising models, and has been, able 
to reproduce all results obtained by other means. The method also sheds 
considerable light on the topological conditions which give rise to an 
exact solution in two dimensions and shows that they are absent in three 
dimensions. However, it is also fair to state that only the knowledge 
of the existence and form of Onsager’s solution enabled Kac and Ward 
to drive their ideas to a satisfactory conclusion. 


tFor further information on this subject see Polya and Szegé (1925), Kac 


et al. (1953), Grenander and Szegé (1958). I am indebted to Professor E. W. 
Montroll for these references. 
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(i) Partition function for a large net. Kac and Ward start from the 
formula (35) giving the partition function for the Ising model in terms of 
closed polygons on the lattice. For a lattice of n horizontal and m vertical 
links connected together on a torus this formula is 


Ziyyl'= 2"(0sh K cosh Kym] 14 > > per. syore |, 


r=1s=1 


(Oe tenh iC wep etanbe | ee, . . 4. (112) 


where p(r,s) is the number of closed polygons with even, vertices which 
can, be constructed on the lattice from r horizontal and s vertical links 
(fig. 7). Following Kac and Ward we shall endeavour to construct a 
determinant each term of whose expansion corresponds uniquely to one 
of the configurations in the sum (112). 

Denote the elements of the determinant by Ap . Then any term in 
det A(= A) is of the form 


+ (Ayp Aop, ... Aypy) 
and can be broken up into cycles 
+ (App, --- Ap,p,)(Ag,o,--- Agya,)--(Arr,--:4z,2,)) + (113) 


We would like each cycle to represent one closed polygon, and we must 
choose the elements Ap, so that only polygons represented in, (112) enter 
in the expansion of A. Let us take for the suffixes P, Q the links in the 
lattice, but directions must also be specified if we are to trace out polygons ; 
there are four links directed towards each point (i, 7) of the lattice the 
directions being right R(—), left L(<), up U(*t ), or down D( ¥ ). Thus 
each suffix P corresponds to a directed bond in the lattice (7, 7, I) where 
T has four possible values R, L, U, D; A is then of order 4mn. We take 
App: to be zero unless P, P’ correspond to nearest neighbour lattice points, 
i.e. unless 
(tet leo Or t=1, 7=9 +1. 


We now notice that the determinant counts oriented polygons, both 
clockwise and anticlockwise orientations being possible, and being counted 
separately. Kac and Ward were impressed by the fact that in the formula 
(101) of Kaufman, for the p.f. (partition function) of a net wrapped on a 
torus m and n do not enter symmetrically and m/2 enters naturally rather 
than m; thus the square of the p.f. might prove considerably simpler than 
the p.f. itself. If we now square Z,,,, in (112) we can represent this as a 
superposition of two independent sets of polygons on one another, and 
we can take one set to be oriented clockwise and the other anticlockwise 
(fig. 22). 

Now in A every directed link occurs exactly once in each term, whereas 
in the polygon expansion only links forming parts of polygons occur. 
This difficulty is readily surmounted by taking App=1, and hence making 
all the unused bonds in the polygon expansion correspond to single cycles 
in A, 
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We specified above that App. is only non zero when P and P’ correspond 
to nearest neighbour lattice points. But we can be more particular now 
in regard to the possible direction I” of P’. Thus if i’=i+1, j=) then if 
must be R (——>); the only non-zero values can be represented graphically 
as 


| eres 


—. .——> °> e—-> 
immediate reversals 
——_ <—__—_ 
Ores . 


must be excluded and the corresponding elements are zero. Hence of the 
16 terms in the 4x4 section corresponding to (i,j) and (i+1,) only 
three terms are non-zero, A(i,j,R; i+1,7,R), A(i,j,D; i+1,j, R), and 
A(i,j,U; i+1,j,R). Analogous results can be derived for the other 
nearest neighbours of P. 


Z mn. 8 represented by a superposition of two sets of polygons the first oriented 
clockwise (full lines) and the second anticlockwise (broken lines). 


The actual values of the elements of A must now be chosen so as to 
give the correct factor v's to the term in A corresponding to p(r,s). 
Clearly any term corresponding to a given link must have the appropriate 
v or v’ (e.g. in the above example all the terms will contain the factor v’). 
But the terms in A are counted with + sign according as the corresponding 
permutation is odd or even, and we may therefore have to insert an addi- 
tional factor to ensure that terms corresponding to closed polygons are 
always counted as positive. Now it is well known that a cycle with an 
even number of terms corresponds to an odd permutation and hence to a 
minus sign in A; thus the factor must be chosen so that an additional 
factor — 1 accompanies every even closed cycle. It will easily be verified 
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that this can be achieved by taking a factor exp (zi /4) for elements involving 
an anti-clockwise turn, of 7/4 


e e ee 
e.g. | 
— ° ° 
and exp (—7i/4) for those involving a clockwise turn, of 7/4 


——> 


: , 
ie 


e.g. 


other elements remain unchanged 


e.g. 


! 

| 

Sie. Pee e ty 
gt 


For example the portion of A corresponding to (i,7) and (i+ 1,7) is 


1 £6, (i 8. 
R vy’ Os O20 
LL 0 0 0 0 ; 
O(1,0= 7] 2. bog | lanexp (ai/4)]. 
D wa 0 0 0 


Similarly the three portions corresponding to (7,7) and (t—1,7), (7,7 +1), 
(i,j —1) are respectively 


0 0 va O 0 0 Or0 
Oar OF Gas 20) ee oie ea 8. 
0(0,1)= 0 0 v rt) ON) De Geode!) 
0 0 0 0 OF var Ok” 0 
On O80" atey—* 
OR Osan, 0o 
Sect 080.0208 0 
OF-O0070 v 


pare sh (114) 


The construction of the appropriate determinant is now complete. For 
graphs consisting of simple closed polygons with no point or side in common, 
it is clear that the counting performed by the determinant is correct. 
But when we begin, to consider graphs containing polygons with points 
and sides in common, substantial difficulties remain in, establishing that 
the counting is always correct. In fact the terms in the expansion of the 
determinant do not correspond identically with those of the oriented 
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graphs, but in all cases which have been investigated in detail the results 
are equivalent. Sometimes the determinant counts additional terms 
which cancel out in pairs (fig. 23a,b); and sometimes it does not count 
the correct polygons but other configurations which are equivalent (fig. 
23c,d). (Several additional examples are given in Newell and Montroll 
1953). But we shall now show that the partition function obtained from 
A is identical with that derived by Onsager, so it seems clear that such 
cancellations and rearrangements always occur exactlyf. 

The matrix A is clearly cyclic in both dimensions all lattice points 
being equivalent. Let us denote any eigenvector of A by X; this can 
be divided into mn sections X; ; each of length 4 units corresponding to the 
lattice point (i, j). The eigenvectors of A are then given by X; ;=ae'” 
where ¢ is an nth root, and 7 an, mth root of unity, and a is a 4 x 1 vector. 


Fig. 23 


There is not always an exact correspondence between terms of A and Z,,,,” 
but a correspondence of final results. (a), (b) represent curves absent 
in Z,,,” Which cancel in pairs in A. Only (c) are present in Z,,,,” 
whereas A counts only (d). 


The matrix A thus reduces to the 4 x 4 matrix 
I,+ €O(1, 0)+¢—O(—1, 0) + 7O(0,1)+-10(0,-1), . (115) 


and its determinant is found by multiplying the mn 4x 4 determinants 
of (115) for n possible values of « and m of 7. We find that 


Ke II II {(1+ v?)(1 + v0’) + 2v’(1 — v?) cos (277/n) 


r=1 s=1 
+ 2v(1—v'2) cos 2s7/m)}, ‘eo eo ON en 
and hence for large m, n 


n m 
Zinn = 2" TT TI] {cosh 2K cosh 2K’ + sinh 2K’ cos (2rz/n) 


r=1 s=1 
+ sinh 2K cos (2s87/m)}. o> OS SA oe ee 
en ep ee a 
tThis has been established rigorously recently by S. Sherman following 


some ideas initially put forward by R. P. Feynman. See “Combinatorial 
Aspects of the Ising Model for Ferromagnetism ”’ to appear in J. math. Phys. 
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Taking logarithms and proceeding to the limit as m, n->0o we re-derive 
for the partition function per spin the Onsager formula 


In (A/2) = mal { In (cosh 2K cosh 2K’ — sinh 2K cos w — sinh 2K’ cos w’) 
oJ 0 


des des. eee este (11s) 
(ii) End effects. We have assumed throughout our work that end 
effects can, be ignored, i.e. that changes in the partition, function, arising 
from finite boundaries become negligible for large m and n. Temperley 
(1957) has shown that this can be established rigorously. The counting 
determinant for a finite lattice is only different from A in lattice points 
corresponding to the boundary, but its cyclic character is destroyed. 
Temperley shows how to combine cyclic solutions so as to satisfy the 
boundary conditions, and concludes that even at the critical point the 
changes in, Z,,, are only or order m or n. 


(ili) Hxact partition function for a torus. The exact partition function for 
a finite torus can, also be established by the combinatorial method, as was 
shown, by Potts and Ward (1955). The error introduced in, the counting 
determinant above for a finite net arises from graphs which loop the torus, 
and which may be counted with the wrong sign. The wrongly counted 
graphs always include one or more of three simple loops 

(a) a horizontal loop with n bonds counted with sign, (— 1)"-1, 

(b) a vertical loop with m bonds counted with sign (—1)”-4, 

(c) a combination of these two with (m+n) bonds counted with sign 

( as UB 2 setae 


It is evident that a restriction to even or odd n, m is not sufficient to correct 
these signs. However, by replacing K and K’ by —K, — K’ the horizontal 
loops will now be counted with sign (=1)"-? (—1)"=-—1; similarly all 
the other loops will be counted with sign—1. Our purpose must now be to 
correct the negative sign, for all such loops. 

Potts and Ward made use of the following device. Let us denote A by 
P,(K, K’). Consider the counting determinant which changes the sign 
of all interactions in the lattice except one row of vertical interactions 
(fig. 24.2); let us denote this by P,(— K, —K’). Then this counts correctly 
all graphs without loops, vertical loops are counted with a positive sign, 
horizontal with a negative sign, and vertical and horizontal loops with a 
positive sign. Define similarly P;(—K,—K’‘) for the counting deter- 
minant which changes the sign of all interactions in the lattice except one 
column of horizontal interactions (fig. 246); and P,(—K, —K’) for the 
corresponding determinant leaving both a row anda column unchanged. 
Then the counting of loops is summarized in table 1. 

At high temperatures we must count all loops with a positive sign. 
But at low temperatures we use the duality transformation (§ 3.4.2) ; 
from eqn. (38) we notice that a factor 2 enters, because of the degeneracy 
of the lowest energy state, but configurations which loop the torus singly 
do not give rise to a low temperature configuration and must not be 
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Hef 


(a ) (b ) 


(a) P,(—K, —K’) corresponds to a change in sign of all vertical interactions 
except those heavily lined ; (b) P;(—K, —K’) corresponds to a change 
in sign of all horizontal interactions except those heavily lined. 


Table 1. Counting of Loops on a Finite Rectangular 
Lattice Wrapped on a Torus 


Twpe of sraph Counting determinant be. 
yp grap Ly tae Fs t df ea dy a fe aed 


Vertical loop — oo - a 0 


Horizontal loop — “ ie 0 


Vertical and horizontal 
loop 


Other graphs 


counted. It will be seen from table 1 that all the counting will be correct 
for the expression. 

Z=h{+\/P,(—-K, —K')+/P,(-K, —K’')++/P,(-K, —K’) 

+4/P,(—-K, —K’)} oi ww, ame ive i at ew 

where the + sign is taken at low temperatures and the — sign at high 
temperaturesf. 

tIn fact for a finite torus the p.f. is continuous and is adequately specified by 
the high temperature expansion. However, it is a useful check to verify the 


formula directly from the low temperature expansion and the duality trans- 
formation. 
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We must now consider the evaluation of the P, j(—K, —K’). We noticed 
that P,(—K, —K’') was cyclic and could be reduced by taking 
Kit, p=eX, ve) X;, p41 = 9%, i] ( e"=7™"= 1). 0 o (120) 


Although P,(—K, — K’) is not cyclic it can be reduced in an analogous 
fashion. We write 


Xi 541 = 0%; iD j#M, 
Xii=—-8X%,m, 8" =—1, 220X121) 
XKi1, j =7X%,;, p 


and it will be found that the matrix A again reduces to a 4x4 matrix 
with 6 replacing «. The basis of the method will become clear from a study 
of the eigenvectors of the matrix 


Qik @ 0 0 

(Ara 0 0 

ee a eee 
0 0 0 @ il 

=i Ove "0; 0 0 


Hence we can write 


P,(-K, —K’')= II IT {cosh 2K cosh 2K’ — sinh 2K’ cos (2r7/7) 


r=1 s= 
—sinh 2K cos [(2s — 1)z/n]}. Sep eae Ase (E22) 
Similarly P3(—K, —K’) retains « and replaces 7 by « where «™_~1; and 
P,(—K, —K’) replaces « by 6 and y by x. 
We can, write P,(—K, —K’) in the form 
P,(—-K, —K’)=2?m II II {cosh 2K cosh 2K’ — sinh 2K’ cos (2rz/n) 
r=1 s=1 


— sinh 2K cos (2sz/m)} 
=(2sinh2K)"" [] 4sinh? (my,,/2) eee et, (23) 
s=1 


where the y, are as defined in (92), i.e. 
cosh y, = coth 2K cosh 2K’ — cosech 2K sinh 2K’ cos (rz/n). 


Here we have made use of the identity 


Vi Il (cosh y — cos 2r7/m) = 4 sinh? (my/2). . . (124) 
r=1 

Proceeding similarly with P,, Ps, P, and defining y,, appropriately it 
can, be shown that Z,,,, a8 given by (119) reduces to (101) as given by 

Kaufman. 
(iv) The degree of order. Potts and Ward (1955) also showed that the 
correlations between, spins on different sites, and hence the corresponding 
degree of order could also be derived readily by the combinatorial method. 
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If # is the Hamiltonian for the system then the correlation between 
spins on sites (2,2) and (a+ 2,b+2)7 is given by 
Zinn F2, 2Fa+2, b+2) = > Oo, 2Fa+2,b42 0XP (— # kT’). . (125) 


lattice 
sites 


The exponential in (125) can be expanded as in (33) and the terms repre- 
sented graphically as in fig. 6. In the absence of a magnetic field the only 
graphs giving non-zero contributions will be closed polygons with even 
vertices as before, except that at the points (2,2) and (a+2, b+2) there 
must be free ends with an odd number of vertices (fig. 25). Potts and 
Ward converted this into a closed polygon, problem of the type considered 
previously by introducmg extra bonds, so as to close the free-ended graph. 


Fig. 25 


The only graphs to be considered for the correlation between P and Q are 
those with an odd number of vertices at P, Q and an even number 
elsewhere. These are converted into closed polygons by the introduction 
of ‘ extra’ bonds (broken lines). 


The extra bonds (shown, dotted in fig. 25) consist of @ horizontal bonds 
from (2,2) to (a+ 2,2) and 6 vertical bonds from (a+ 2, 2) to (a+ 2,b+42). 
We may now write 

Zinn (2, 2% a+2, b+2) = (2 cosh K’ cosh K)"”" coth® K’ coth? K > x(r, s)v'"v’, 

(126) 

where z(r,8) is the number of closed graphs as before, but counting only 
graphs which include the extra bonds. 

The counting is achieved as before by superimposing the lattice on 
itself, and taking one set of polygons as oriented clockwise and the other 


tThis choice is for convenience to avoid boundary effects in the diagrams. 
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as oriented anticlockwise. This means that the extra bonds must be 
taken in both directions. The extra bonds are taken into account by 
bordering A with 2(a+b) rows and columns corresponding to the 
two possible directions for each extra bond. The elements are chosen 
very much as before, but the diagonal elements of the bordering section 
are taken to be zero, so that these extra bonds must be included in any 
non-zero term. 

The details of the construction and reduction of the determinant are 
given by Potts and Ward. The complete identification of their results 
with those of Kaufman and Onsager, appears to be somewhat involved, 
but calculation in particular cases showed agreement. The method certainly 
provides the simplest practical approach to the evaluation of correlations. 


Fig. 26 
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VA 
A SVAVAVA 
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The triangular lattice is topologically equivalent to a rectangular lattice with 
one second neighbour interaction. 


3.5.4. Exact solutions for other two-dimensional lattices 


(i) The triangular lattice. We have seen, in §3.4.3 how the duality and 
star-triangle transformations enable the Curie point of the triangular 
lattice to be located. Shortly after the publication of the work of Onsager 
and Kaufman, several papers appeared independently giving the exact 
solution for this lattice (Houtappel 1950, Husimi and Syozi 1950, Newell 
1950, Temperley 1950, Wannier 1950). Potts (1955) applied the combina- 
torial method to the triangular lattice, and this leads most directly to 
thesolution. From this point of view the triangular lattice is best considered 
as a rectangular lattice with a single second neighbour interaction (fig. 26) ; 
the lattice can be wrapped on a torus and the procedure is then closely 
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analogous to that of §3.5.3. The resulting partition function per spin is 
given (for general interactions K, K’, K") by 


Qa (20 
In (A/2) = (822) | | In {coro + 88'S" —Scosw, 
0/0 


— S’ cos w,—S" cos (w, + w9)} dw, dws (123) 
C;0',C" =cosh2K, cosh2K’, cosh2K”. 
S,S', 8S” =sinh2K, sinh2K’, sinh2K”. 


It is clear that as K”->0 we regain the formula (94) for the rectangular 
lattice. 

The physical properties of the ferromagnetic triangular lattice are 
very similar to those of the quadratic lattice. The singularities in the 
thermodynamic functions are of the same type, and when K=K’=K" 
these functions differ only slightly from one another in relative value. 
As an example of the order of magnitude of this variation the approximate 
formulae for the specific heats in the neighbourhood of the Curie point for 
the two lattices are as follows: 
~ —0-4991In |7’—T,| (triangular), \ (128) 


GIN é 
~ —0-4945In|7—T| (quadratic). 


(Houtappel 1950). 
An interesting and novel feature of the triangular lattice is the anti- 


ferromagnetic solution, and this aspect is discussed in detail by Wannier 
(1950). Antiferromagnetism does not fit into a triangular lattice and we 
can readily show that the energy of the lowest state is only one-third 
of the corresponding ferromagnetic value. First we note that the energy 
cannot be lower than this value ; the net can be completely decomposed into 
triangles each of which contains three interactions. The best we can 
achieve within each triangle is to have two spins of one kind and one of 
another; hence the minimum possible energy is VJ(J <0). 

We now construct states with exactly this minimum energy. Figure 27 
which consists of alternating rows of positive and negative spins is an 
example of such a state. Figure 28 shows another arrangement in which 
the rows can be fitted together at random without changing the energy ; 
the weight of states of this kind is of order NV", and they do not, therefore, 
give rise to any entropy. In fig. 29 we have an arrangement with a finite 
entropy. In this arrangement the correct energy is obtained by an un- 
balance ; each + spin is surrounded by 6—spins, but each — spin contains 
3 + spins and 3 — spins among its neighbours. From this it follows that 
each of the encircled negative spins can be overturned independently with- 
out changing the energy; there are thus at least 2/* independent states. 

But this does not exhaust all possible states. It may happen that a 
trio of neighbouring encircled spins have the same sign (in fact each 
possible configuration of these spins will occur on the average on one-eighth 
of the occasions) and that the spin at the centre of the tripod then has 
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equal numbers of + and — spins as its neighbours, and can be reversed 
without changing the energy. 

The number of states is hence increased by a factor 24/2, Thus the 
entropy at absolute zero is greater than 


2 Nkin2=0-288811Nk.. . . . . . (129) 


Fig. 27 
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A minimum energy anti- ferromagnetic state for the triangular lattice. Rows 
of positive spins alternate with rows of negative spins. 


Fig. 28 
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An arrangement of minimum energy in which rows can be stacked at random. 


Wannier derived an exact expression for this entropy in the form 
7/6 
S(0)/Nk= =| In (2 cosw) dw =0-32306. . . . (130) 
Ti Jic0 


He also pointed out that long-range order was incompatible with a 
degenerate lowest energy state of the kind described. As seen, in fig. 30 
it is possible for two states with different order to adjoin being connected 
by a transition region with the same energy as the ordered regions. It is 


P.M.S. Q 


214 C. Domb on 


not surprising therefore that the partition function for an antiferromagnetic 
triangular lattice (in the absence of a magnetic field) has no singularity, 
and the specific heat curve is smooth and continuous with no infinity. 


Fig. 29 
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An arrangement of minimum energy having finite entropy. The encircled 
spins can point in either direction. 


Fig. 30 
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Two adjoining regions having different types of long range order. The energy 
of the transition layer may be made just as low as that of the ordered 
regions. 


Since the presence of a small magnetic field at 7’=0 would cause a large 
change in the number of spins pointing in the direction of the field, we 
should expect the susceptibility of the triangular lattice to be infinite at 
T’=0. Insome respects, therefore, the behaviour of the antiferromagnetic 
triangular lattice in zero field is analogous to the behaviour of an antiferro- 
magnetic linear chain, in a critical field (§ 3.2.1). 
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(ii) The honeycomb lattice. The solution for the honeycomb lattice can 
be obtained by applying the star-triangle transformation to the solution 
for the triangular lattice. The partition function is given by 


Qar (2a 
In (A/2) = (1622)-2 | I In }{CO’C" +1—8'8" cose, 
0 0 


—8"S cos w,— SS’ cos (w1—w.)}dw day, . . (131) 


Fig. 31 


fe} 0.2 0.4 0.6 0.8 1.0 1-2 1-4 1.6 1.8 2.0 
ee 
A comparison of the specific heats of the honeycomb and quadratic lattices. 
honeycomb ; ——--—--—--— quadratic. 


The physical properties of the honeycomb lattice are again closely similar 
to those of the quadratic lattice. The specific heat approximation near 
the Curie point is (Houtappe 1950). 

Grin 0476) nll Pele. ei jn es 4 (182) 
The specific heats of the honeycomb and quadratic lattices are compared 
in fig. 31. 

Several properties of the quadratic, triangular and honeycomb lattices 
can, be combined into one general formula. Thus Onsager noted (Wannier 
1945) that the Curie points of all three lattices are given by 

Ca ime (tee he Gn be (183) 
where gd denotes the gudermanian, [gd «=sec~'(cosh a) = tan (sinha) | 
and q is the coordination number of the lattice. Similarly we have seen. 
that the inversion, transformation, can be put in the form (59) which applies 
to all three lattices. Finally Potts (1952) showed that the spontaneous 


Q2 
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magnetization of the triangular lattice could be put in the form (111), the 
appropriate transformation variables being given by the inversion 
transformation (58) of §3.4.3 (K+>K*). We may thus write for both 
quadratic and triangular lattices 
(I/a)8 = 1 —(v- — v)?/ (21 —2)?, Fo As ey eee 
using the v variables of § 3.4.2, § 3.4.3; from (60) the spontaneous magneti- 
zation of the triangular lattice is given by 
1626 

(1 + 32?)(1 —2?)8 

In order to derive the spontaneous magnetization of the honeycomb 
lattice Naya (1954) introduced the idea of a ‘semi-ferromagnetic’ honey- 
comb lattice in which one of the two triangular sub-lattices into which the 
honeycomb lattice can be divided (fig. 11) consists of spins with normal 
magnetic moment », and the other consists of spins with zero magnetic 
moment. Naya showed that the spontaneous magnetization of the semi- 
ferromagnetic honeycomb lattice was exactly half that of the normal 
ferromagnetic honeycomb lattice. Now the non-magnetic spins could 
be removed by the star-triangle transformation, (§ 3.4.3) in the presence of 
a magnetic field (§3.4.5), and a relation could thus be derived between 
the spontaneous magnetization, of the triangular lattice and the semi- 
ferromagnetic honeycomb lattice (see also Fisher 1959a). It was found 
that the spontaneous magnetization, of the honeycomb lattice could also 
be put into the form (134), the appropriate transformation, variable being 
given by (62). Hence for this lattice 


(L/w)8 =1— (135) 


1623(1 +23) 
(L-=2)8(L aa"? 

The spontaneous magnetization, varies only slightly from one lattice 
to another among the three lattices considered so far, as will be seen from 
fig. 32. 

The honeycomb lattice has a Curie point for all values and signs of the 
interaction parameters. 

(iii) The kagomé lattice. In §3.4.4 we showed how the decoration 
transformation could be used to derive a relation between, the honeycomb 
and kagomé lattices in the presence of a magnetic field. Hence the thermo- 
dynamic properties of the kagomé lattice can be derived from those of the 
honeycomb lattice (Syozi 1951, Naya 1954). The kagomé lattice does 
not satisfy the relations (59), (133), (134) derived for the other two- 
dimensional lattices. Its antiferromagnetic properties are similar to 
those of the triangular lattice (Kano and Naya 1953), and it possesses a 
finite entropy, 0-50183Nk at the absolute zero. 

(iv) General treatment of two-dimensional lattices. A general approach 
to a large variety of two-dimensional lattices, which includes all the lattices 
discussed so far, and can be applied to many others, was put forward by 
Utiyama in 1951 (see also Syozi 1955a). The methods of Kaufman and 


(I/2)8=1— (136) 
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Onsager can equally well be used for a chequer lattice formed by replacing 
all the black squares of the chessboard by one and the same elementary 
figure. The elementary figure consists of a rectangle with v vertical 
crosspieces (fig. 33). By varying the interaction parameters, and taking 


Fig. 32 
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The spontaneous magnetizations of different lattices. —-—-—----—- honeycomb; 
quadratic’; -++ +--+: triangular. 


Fig. 33 


Elementary figure suggested by Utiyama to be placed on all black squares of 
the chessboard. 


some as 0 or 00, we can generate different lattices. For example with 
y=0, J=oo, J,=0, or J=Jy, Jo=Jo' gives the quadratic lattice; J=0 
gives the honeycomb lattice; J = o gives the triangular lattice (fig. 34). 
With v=1 J= oo corresponds to the kagomé lattice (fig. 35); and Jg=0 
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gives the 4-8 lattice (fig. 36) which is of interest since, although all points 
are equivalent, all bonds are not equivalent ; in this respect the 4-8 lattice 
resembles the three-dimensional close-packed hexagonal lattice. 

The partition function per spin can be written down quite generally 
in terms of the interaction parameters. Thus for y=0 we can write, using 


a notation analogous to (127), 


l 1 20 (27 . : : c 
nWV2)= 55 | | ~ W2[OCyC)Oy + S84 SpSy/ +1 


ss Bodo. cos (w, + Ws) a SS, cos (w, — Wy) —_ (SSp + SS) cos Wy 


Fig. 34 


The chequer lattice with »v=0 gives rise to the quadratic (J =J,, J>=J)’), 
honeycomb (J =0), and triangular (J = co) lattices. 


where n,, is the number of bond strengths J, J,, Jo, Jo’ which become 
infinite, the corresponding linked spins being regarded as one. A similar 
formula was given by Utiyama for general v. As a particular example, 
the partition function for the kagomé lattice is given by 


22 (20 
In (A/2) = (24m) | | In J{O® + S84 20383 + 302 — 2(C'S3 + C282) 
0 0 


x [COS wy + COS Wy + COS (w, + We) ]} dw,das, Se ae. 
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Fig. 35 


Jo Ja fein atk 
The chequer lattice with v=1 gives rise to the kagomé lattice (J = oo) and the 
4-8 lattice (J,=0). 
Fig. 36 


The 4-8 lattice. All vertices are equivalent but bonds belong to two groups, 
square and octagonal. 
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The approach of Utiyama indicated that the general formulae (59), 
(133), (134) apply only to lattices corresponding to v=0; they break down 
for the kagomé lattice for which yp=1. The formula (133) of Onsager which 
determines the Curie point was also generalized to all lattices with v=0 


in the form 
gd2K+9d2K,+gd2K,+gd2K,=7.. . . . (139) 


Using the above approach and the duality, decoration, and star-triangle 
transformations, Syozi (1955b) determined the Curie points of many 
complex lattices, some of them having average coordination number 5. 
By increasing the number of decorating spins it was possible to loosen 
the coordination of the lattice until the Curie point ultimately tended to 
an); 

We may conclude generally that the mathematical nature of the 
singularities is the same for all of the two-dimensional lattices considered, 
and that the detailed thermodynamic properties of the model are very 
insensitive to change in lattice type. 


3.6. Hxact Series Expansions 


All the results which we have derived so far correspond to two-dimen- 
sional lattices in the absence of a magnetic field. By applying certain 
transformations to the Onsager solution Fisher (1960) has derived an 
exact solution for an antiferromagnetic in a magnetic field; this corres- 
ponds to a specialized model, and we will defer consideration, of its detailed 
properties to a subsequent publication. The methods described so far 
in this chapter are inapplicable in the presence of a magnetic field or in 
three dimensions. The most important source of exact information regard - 
ing the latter problems are series expansions of the partition function at 
low and high temperatures. 

The derivation of substantial numbers of terms of these series gives 
rise to topological questions which will be considered in detail in §5. 
Our present purpose is to formulate the problem, and establish the form 
of the expansions in different regions of the (u,z) plane. We shall also 
find that an investigation of the propagation of order, and of the presence 
or absence of long-range order, will help us to identify the singularities 
of the partition function. 


3.6.1. Low temperature expansions 


We have already noted (§ 3.4.2) that at low temperatures the partition 
function can be expanded in terms of successive deviations (corresponding 
to more and more overturned spins) from the lowest energy state. Let us 
consider a set of V interacting spins (J > 0) in a magnetic field in direction 2. 
(It is more useful in this case to derive a low density expansion.) Following 
the discussion and notation of § 2.2, the lowest energy state consists of all 
spins pointing in direction 2; there are YO, excited states corresponding to 
one overturned spin, *C, corresponding to two overturned spins, . . . ¥C, 
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corresponding to r overturned spins. The “C,, spins can be divided into 
different energy groups; the highest energy deviation from the ground 
state, (2rmH + 2qrJ), arises from r spins, no two of which are nearest 
neighbours; the next from r spins with one nearest neighbour link, the 
next from r spins with two nearest neighbour links, and so on. Each 
nearest neighbour link decreases the energy by 2; as an example the 
configurations which correspond to r=3 are shown diagrammatically in 
fig. 37. We may thus expand A,(,z) in the form 


Ay(p,2)=1+pF\(N,2)+p?F(N,z)+...p"F(N,2) + ---(u<l) (140) 


where there are “C, terms contributing to F,(N,z). For a simple quadratic 
lattice wrapped on a torus it is easy to see that 


N(N-5 

PN, 2)=N2*, BN, 2) =2N28+ ee 

141) 

Bie 2 P) ( 

F(N,2)=6Nz + 2N(N — 8)z10 + eee 

Fig. 37 
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Low temperature configurations with three spins. 


So far the expansion has been exact. We shall now allow N to become 
large and derive a corresponding expansion for A(u,z). We shall first 
deduce certain relations formally, and will subsequently discuss the 
detailed justification. For large N we should expect on physical grounds 
that 

DAMON UIEN ean ieee toa ce es 1 (142) 


there might be a multiplying factor whose effect would become negligible 
for large N. Hence if we assume that A(,z) can be expanded as a power 
series in x of the form 


A(u,2)=14 wfi(2) + fez) +-.-wh(Z) +... (u< I) . (148) 


we can equate successive powers of NV, and deduce the values of the /,(z) 
from the F,(N,z). In fact, for a lattice in which the edges are connected 
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cyclically, there are no surface effects, and it will be found that the multi- 
plying factor is unity. We then have an identity 


1+ pF (N,z)+ pw? P(N, 2) +. Us Pol Ne) oe 
=14 °C, [nAle) +e? fol) +.-- whl) + 
+O, [uf (2) + wfelz) + --- F(z) + 


C [pf (2) 2s. Se f(a) eae i NYAS} 
Putting NV = 1, we find the very simple relation . 
f,{2) =F, (2) Se eee te Oe Ee ee ee 


(Hence for the 8.Q. lattice 
filz)=24, fa(z) = 22®—228, f(z) = 628 — 14279 + 82'). 


The derivation, of this result depends on, the identity of both sides of 
(144) and to convince oneself of its correctness it is advisable to write out 
fully the first few terms of each side of (144) in order to verify that this is 
an, identity. 

It is sometimes more convenient to use the expansion for In A(p,z) in 
the form 

InA(y, 2) =pgilz) + w2gol2)+---wg(2)t---. (<1)... (146) 
the g,(z) being closely related to the cluster integrals in the Mayer theory 
(see § 4.5 for further details). If we take the exponential of (146) we 
can, relate g,(z) to F,.(N,z); expanding F’,(.V,z) in powers of V, we can 
write 

EN, 2) = Ayy(Z)N 6+ Gip3)(Z) NT +... eg (ZN. 2. (147) 


We now formally expand exp [NV In A(u,z)] in powers of the variable NV 
and equate coefficients of powers of NV; we deduce that 


g,(2)=a,4(2)i kos 12.) oot See ee 


that is g,(z) is the coefficient of NV in the expansion of F',(V,z) in powers 
of VN. (Thus for the 8.Q. lattice 
gi(z)=24, go(z)=225— $28 ga(z) = 628 — 1621 + 104212), 

The justification for the results obtained above lies in the consideration 
of the form of F’,(V,z) for increasing values of r/N. When the concen- 
tration of overturned spins is small it will be rare to find any doublets 
(pairs of adjacent overturned spins), and a treatment which pea 
only singlets will provide a reasonable approximation. Choosing /;(z) 
satisfy (145), and ignoring all higher terms of the series (143), pate we 
to taking 

FN, 2) =O, f,(z) d+ Chea, Hee 
and hence to neglecting doublets in FN, z). eae stopping the series 
after two terms corresponds to taking account correctly of doublet forma- 
tion in F’.(N,z) but not of higher clusters. Hence we see that the higher 


terms become significant as r/N increases, and can be dealt with in the 
same manner. 
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The argument is similar to that used by Mayer in the derivation, of his 
cluster integral expansion (Mayer and Mayer 1940). For large N the 
series for Ay(y,z) can be replaced by its maximum term; the evaluation 
of this maximum term as a function of ~ and z would lead to relations 
equivalent to (142), (145) and (148). 

We have so far ignored states in which more than half of the spins of 
the assembly point in direction 1, and in the presence of a magnetic field, 
however small, this is certainly valid for large N. When the field is zero 
these states become significant ; we can easily take account of them because 
of the symmetry between directions 1 and 2, and we can write instead of 
(140) 

Ay(ps2) = (1+ pe) + (ut pA), (N, 2) + (u2 + pk 2)P(N, 2) +... 
Bee tepie ws EEG) tomers (forsadl y1)e tse. 2. (150) 


It will be verified that the relation (28) is now satisfied by Ay(u,z). When 
pis less than 1 the second term in each bracket of (150) is negligible and 
we recover (140). 

When p=1, both terms in each bracket are equal, and the expression 
(140) must be multiplied by a factor 2. On taking the Nth root, the 
expression (143) for A(u,z) must be multiplied by 2', and this differs 
negligibly from unity for large VN. Hence A(y,z) is continuous at w=1. 

The magnetization is given, by 


0 0 
Ly = kT pln Zy= | N24 Im Anly2) | ee ae) 


When » = 1, it will be seen from (150) that 0/.y/¢u is equal to NA y, and hence 
that J=0. When p<1, however, if we use the expansion (146), the 
magnetization is equal to 
— Nok —2[pgi(2) + 2y792(2) +... rp"g,(2)+..-]}, . «>. (152) 

and, as long as the expansion, exists at .=1 for sufficiently small z, this 
tends to anon-zero value as y>1. Hence there is a spontaneous magnetiza- 
tion for sufficiently small z. The existence of a spontaneous magnetization 
therefore depends on the possibility of a valid expansion, of the type (146) 
for »=1 near z=0. For a one-dimensional chain every f,(z) starts with 
z2, since any number of spins can, be overturned with only two adjacent 
1-2 links; hence the coefficients of z? in the expansions (143), (146) become 
infinite when, »=1 and there is no spontaneous magnetization. Similarly 
for a strip of any finite width, since a finite number of overturned spins 
can divide the chain into two parts, one of the coefficients in the expansion 
would become infinite. It is only when the model extends to infinity in 
two or more dimensions that high powers of » necessitate correspondingly 
high powers of z, the expansions (143), (146) are valid when » = 1, and there 
is a spontaneous magnetization. 

The above discussion, is not put forward rigorously but to enable one to 
understand physically how the discontinuities in magnetization arise for 
the Ising model. j 
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In the antiferromagnetic case (J <0) the expansions (143), (146) are 
still valid for fields sufficiently large compared with the coupling between 
spins. For small fields the same method can be used to derive a low 
temperature expansion for lattices which can be decomposed into two 
equivalent sublattices (§ 2.2). When considering the configurations which 
correspond to overturned spins (fig. 37) we must now differentiate between 
the two sublattices, and the problem is therefore more complicated. By 
analogy with § 2 (9) we write the partition function in the form 

Zyl =y PAANOp,y), [y=e*=exp(2J/kT)] . . (153) 
and A™ (u,y) can be expanded as a power series in y. For example, the 
first few terms of this expansion for the S.Q. lattice are (Brooks and Domb 
1951) 

NOp, y) = 1+ 3(e + w)y* + 2y° 
+ [—(u* + p7)/8+3(u> +p) — Fly*+ .... (154) 
There are no discontinuities in A@(y, y) as 1>1, and there is no spontaneous 
magnetization. As might be expected physically the singularities in an 
antiferromagnetic occur in the neighbourhood of the critical field, and are 
caused by the transition from the state of alternating order, represented 
by (154), to the state of lined-up spins represented by (143). 


3.6.2. High temperature expansions 

We have already noted in § 3.4.1 that eqns. (31), (32) and (33) provide a 
basis for expansions of the partition function in powers of tanh K, i.e. 
high temperature expansions. Expansions of this kind were used by 
Kirkwood (1938) in considering order—disorder transitions, and by 
Opechowski (1937) for the Heisenberg model; they were further developed 
by van der Waerden (1941), Kramers and Wannier (1941), Oguchi (1949), 
Wakefield (1951), Sykes (1956) and others. 

The detailed discussion is closely analogous to that of the previous 
section with line configurations replacing spin configurations, and 
replacing N (fig. 6). We must sum every configuration over the o; variables ; 
any one of the o; which does not occur explicitly in a term of the expansion 
(33) gives rise to a factor (u-"?+ 2") in this term; likewise any o; which 
occurs to an even power. Any o; which occurs to an odd power gives rise 
to a factor (u-!*—."?). Thus we may write 

Zyi = (u-V? + w?)¥ (cosh K) [1 + tanh K O,(4,r) + 
+tanh? kK ®(/%,7)+...tanh"K ®,(W,7)+...], vats SLORY 
7=(1—p)/(1 +1). 


Putting /=qN/2, and comparing with (140), we derive an expansion 
for the partition function per spin in the form 


Z= (wu? + w'?)(cosh K)4?[1 + tanh K ¢,(7) + tanh? K g,(7) 
+...tanh’ K ¢,(7)+...] tial> tists ie ve ae COO 
where, following (145) 
Pelt) = O,(g/2, 7). LS ete (EBT) 
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This expansion is valid for all ., and has smooth analytic behaviour at 
#=1, the magnetization being zero on this line. It is also applicable to 
antiferromagnetism (K < 0). 

High temperature expansions can, also be thought of as arising from the 
moments of g(V; N,,N,,) regarded as a probability distribution. We 
have from § 2 (5) and (9) 


Ay(u; 2) = > g(N; Ny, N yo) pie 


Ni, Nis 
ice) T s(9R\8 
where ZNu — exp (—2N,,K)= > (-y (158) 
« t(9 7, \t 
wt =exp(—2N,L)= > (—1) a (L= mH |kT) 
t=0 : 
Hence we can write 
2 N,®N,') (2K )8(2L)¢ 
NAA; z)= Py (— 1)s+ < 12 ae ) ( ) ‘ (159) 


and Ax(y,2) is the moment generating function for a completely random 
mixture. The properties of a random mixture for a simple quadratic 
lattice are of statistical interest in relation to the infection of trees in an 
orchard and were considered by Moran (1947). He calculated several 
moments, and proved that the distribution is normal for large VN. The 
latter result is not surprising physically, and tells us that 


In Ay(u, 2) ~ NV In A(p,2) 


for large NV; for detailed knowledge of the partition function and its 
singularities we require exact values of all the cumulants of the distribution. 


3.6.3. Fluctuations 
The expansion discussed in §3.6.1 divides the partition function into 
terms corresponding to various possible configurations of the assembly, 
and by Boltzmann’s principle each term is proportional to the probability 
of achieving the corresponding configuration; thus from (150) the 
probability of exactly r spins pointing in direction 1 is p’P,(N,z)/Ay(u, 2). 
Hence the mean number of spins pointing in direction | is given by 
My )Ay (#2) = Ne + [1+ (ND PN, 2) 
+... [rp +(N —r)p* P(N, 2)+.... a) 


When, 1. < 1 the second term in each factor on the right-hand side can be 
neglected, and 
Ay (Ny » yes OAy/Op = NpA*-1 OA/Op. 
CONG aN wit a. eo Sera (h61) 
However when p=1, (N,)=N/2, and there is a discontinuity in (N,) 
at p= 1. 
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Similarly (N°) is determined by (u0/0u)’Ay(u,2). When p<1 the 
mean, square deviation is given by 
(N,2)— (NP = [pA JA t+ pA, /A—(wA,/A)*], - (162) 
and the fluctuations are normal. But when »=1 we obtain instead 


7 


(Ny?) — (Ny)? = a [1 — 20(z) }? + N[8(z) + A(z) —6(z)"], . (168) 
where 6(z)= Lt wA,/A, 8(z)= Lt p?A,,/A. Therefore, as long as @ is 
ul 


ul 
less than half the fluctuations in V, are greater than, normal, the r.m.s. 
deviations being of order NV. Fluctuations in V,, can be determined by 
differentiating partially with respect to z instead of «4; the results can 
also be confirmed from eqn. (159). 


3.6.4. Haupansions for the degree of order 


The methods of the previous sections can be used to provide expansions 
at high and low temperatures of the degree of order discussed in §3.3 
(Oguchi 1951). At high temperatures we have already defined the appro- 
priate configurational problem for the quadratic net in (125), and this 
can readily be extended to other cases. We have the corresponding general 
relation for the correlation beween spins a and b 


oe) i 


the detailed expansion of which can, be undertaken in the same manner as 
(155). The first term in the expansion, arising from no lines, will give rise 
to a term 7”; this is equal to (c,,) oc, ) to the first order, and corresponds to 
the absence of correlation, at sufficiently high temperatures. Correlation 
is introduced by line configurations (e.g. fig. 6), and ce”, can be determined 
from (19). 

At low temperatures, for a ferromagnetic in a magnetic field, we can 
determine (oa,,;“") (§3.3) by a method similar to that of §3.6.3. Let 
F,.“(N,z) represent the sub-group of F,.(.NV,z) in which spins a and 6 are 
both overturned. Then 

Ay(H4,2)<ay,) = p2F (WN, 2) + pF (N, 2) +... pF O(N, 2) +... 
(165) 

{a,) is equal to (N,)/N and is therefore determined by (161). Hence we 
can calculate c/,,, from (18). In the presence of a magnetic field there is 
only one lowest energy state; correlation between the specified spins 
arises only from configurations with continuous chains of overturned 
spins containing both specified spins. If the two specified spins are 
sufficiently far apart no such configurations enter, (a4;°) =(a,)2and o% , =0. 

In, the absence of a magnetic field, however, both of the equivalent lowest 
energy states must be taken into account. If the lowest energy state at 
absolute zero consists of all spins pointing in direction 2, the mean fraction 
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of spins pointing in direction 1 at any temperature is equal to 6(z); from 
(161) and (163); the mean fraction pointing in direction 2 is [1 — O(z)]. 
Thus since both lowest energy states are equally probable, 


Cor” )= 3 [0(z) P+ [1 — O(2) P= 3 [1 — 20(z) + 20(z)?]. (166) 
{a ) is equal to 4, and hence from (18) 
(on) =H te? 4) =4 (1?) 167) 


_ Thus we deduce that 

Torte 00 ecm Wenn <n 2°41 68) 
i.e. that the long-range order is equal to the spontaneous magnetization. 
We see also that long-range order arises physically as a result of the inter- 
ference of equivalent lowest energy configurations. 

For an antiferromagnetic in the absence of a magnetic field the result 
(168) remains valid, although there is now no spontaneous magnetization. 
But in the presence of a magnetic field a new calculation is required since 
long-range order is present for all magnetic fields of magnitude less than a 
critical value. We are concerned with lattices which can be sub-divided 
into two equivalent sub-lattices, so that the lowest energy antiferromag- 
netic state corresponds to the spins of the sub-lattices all pointing in direc- 
tions 1 and 2 respectively. These sub-lattices can be interchanged without 
affecting the energy, and the lowest energy level is therefore degenerate ; 
this last result is true in the presence of a magnetic field, provided that the 
field is not sufficiently large to break the coupling between the spins and line 
them up (§ 2.2, § 2, fig. 2). Suppose that the field points in direction 2, and 
at any temperature 7 let (7,) be the mean fraction of spins of sub-lattice 2 
which are turned over from direction 2 to direction 1, and (y,) the mean 
fraction of spins of sub-lattice 1 which have turned over from direction 1 
to direction 2. Then 

Cay sl) — ep) «(1 89) 


To evaluate (a,,°) we must, as before, take equal account of both lowest 
energy states, and we find 


{017 =3(1— (me)? + 34m) ee ee (170) 
Hence from (18) we deduce that 
BS 1— (m1) — (2) Roy NEO 2<m1) ; eNe(LTL) 
1+ <1) — (Ne) 1+ (m1) — (2) 


Tn the absence of a magnetic field both (7,) and (72) are equal to @(z), and 
we recover (168). It is interesting to note in the general case that R be- 
comes zero when <7, )+ <7.) =1, i.e. when on the average half of the spins 
in the lattice have overturned. 


3.7. Location of Singularities 


Having derived various types of series expansion, and specified the 
conditions of occurrence of long-range order, we are now in a position to 
outline the singularities of the partition function in the (4, z) plane. These 
singularities are associated with the disappearance of long-range order. 
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Quite generally if we have established the existence of long-range order in an 
assembly for a finite temperature range above 7’'=0, we know that a 
singularity must occur at a temperature 7’. at which the long-range order 
becomes zero (for when the long-range order is non-zero it can be used to 
specify the temperature, and the assembly possesses an, additional thermo- 
dynamic variable). For one-dimensional assemblies the long-range order 
which exists at 7'=0 disappears at a finite temperature however small. 


Fig. 38 
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The singularities of the Ising model of a ferromagnetic. All singularities lie 
on AC. In the region «<1 there is no long-range order. 


Let us first consider a ferromagnetic (Domb 1949b). We need only deal 
with the range 0 <p, 2< 1, and this is illustrated in fig. 38. We know that 
in the range y< 1 there is no long-range order, and in fact long-range order 
exists only when »=1, and at temperatures below the Curie point repre- 
sented by C. Thus the only ordered region in the plane is the line AC. 
There is a discontinuity in A, as p> 1 for z< ze. 

The low density expansion represented by (143) corresponds to starting 
at some point of DE and moving parallel to DA. The low temperature 
expansion when y= 1 corresponds to starting at A and moving along AC. 
The high temperature expansion (156) corresponds to starting at some point 
of BE and moving parallel to BA. From our discussion we should expect 
the partition function to have no singularities except on the line u.=1; 
this was established rigorously by Yang and Lee (1952), who considered 
very generally the conditions which give rise to singularities in cooperative 
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assemblies. In the presence of a finite magnetic field, therefore, the 
thermodynamic functions behave smoothly and analytically. We shall 
also find that the two expansions (143) and (156), which we derived by 
independent methods, are closely related and can be deduced from one 
another. 

The antiferromagnetic corresponds to the range 1<z< but we shall 
prefer to use the variable y=z~", so that again 0 < HL, ¥<1 (fig. 39). Fora 
lattice which can be sub-divided into equivalent sub-lattices the solution 
when 4 = 1, represented by A’B’, is identical with that for a ferromagnetic, 
the Curie point being at C’. There is now no spontaneous magnetization 
as p> 1, and the partition function behaves smoothly and analytically in 
this region. However, the long-range order which is present along A’C 


Fig. 39 
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y 
The singularities of the Ising model of an antiferromagnetic. The line of 
singularities separates the region with long-range order from that with 
no long-range order. 


persists in the presence of a magnetic field, and the points in the (, y) plane 
at which it disappears trace out a curve of singularities D’F’C’ dividing 
ordered from disordered states. At low temperatures we should expect this 
curve to correspond to the critical field (§ 2.2), and hence to be given by 
u=K,y* where K, is a constant. In the disordered region the expansion 
(156) remains valid starting from a point on B’E’ and moving parallel to 
B’C’, and the expansion (143) starting at a point on D’E’ and moving parallel 
to D’A’. In the ordered region the expansion (154) can be used starting 
at points on A’D’ or A’C’. There are so far no general theorems for an 
antiferromagnetic of the type enunciated by Yang and Lee for a ferro- 
magnetic; however, it seems reasonable to assume that the demarcation 
of singularities given here is valid. 


P.M.S. R 
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For lattices such as the triangular there is no ordered state in the absence 
ofa magnetic field. However, itis easy to see that when the ratio of positive 
to negative spins is 2 to 1 an ordered state can be constructed, and there 
will be an appropriate critical field for this ordered state. Hence it seems 
very probable that the singularities are as indicated in fig. 40. 


Fig. 40 
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Singularities of the antiferromagnetic triangular lattice. In zero field there is 
no long-range order. In the presence of a magnetic field long-range 
order can exist. 


3.7.1. Application to regular solutions 


We have seen from § 2.3 that the thermodynamic properties of regular 
solutions are completely determined by the function A (u,z). Weshall now 
show how to interpret the singularities for the theory of regular solutions. 
From § 2.3 (10)—(13) the free energy for a solution of concentration o, 
is given by 
FS = —kT [In A- oy In pe] t qayey1/2+9(1— oy) eon)/2  . . . (172) 
where 

MAGIA iy. low. hae eae ene 


The relation (173) determines a path in the (u, 2) plane along which A must 
be evaluated to determine FS. When a,= 4 the path is AB (fig. 38). When 
a is different from } a different type of path results. For z< ze the general 
behaviour of wA,/A as a function of p is similar to the curve in fig. 41. Its 
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value increases steadily from zero at ~»=0 to 6(z) at w=1 and there is a 
jumpto}at~=1. This discontinuity arises for an infinite lattice, and would 
be a very rapid continuous change for a large finite lattice. Hence if 
O(z) < a,< % the relation (173) can only be satisfied on the line AB. For 
z=1 we can use the high temperature expansion (156) to show that 
v/(1+p)>o. Thus for any value of «, different from 4 the path given, by 
(173) moves along AC until 0(z) = «,, and then it turns into the plane cutting 
BE in a point ~=a,/(1—«,). Typical paths are AC,B,, AC,B,. 

Clearly a transition takes place at points such as C,, and we shall interpret 
the portion AC, as corresponding to a separation, into two phases, and the 
portion C,B, as solution. This interpretation is borne out by the curve of 
free energy as a function, of concentration shown in fig. 42. (For simplicity 
the case €,;=€. has been taken, but the resulting features are essentially 


Fig. 41 
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Discontinuous behaviour of A, at »=1 for 7’<T,. The function pA, increases 
steadily from 0 at ~=0 to 6(z) as w 1, then jumps discontinuously to 
ai 
> at = le F 


the same for the general case.) The portion AC, may be regarded as a 
mixture of two phases at Y, phase I being a solution of 6(z) of substance 1 
in substance 2, and phase II a solution of 6(z) of substance 2 in substance 1. 
6(z) is thus the maximum concentration which can exist in solution at the 
given temperature, and is in fact the solubility curve for the two substances. 
As we move along AC, the proportion, of phase II gets smaller until at C, it 
disappears. The nature of the transition at Cj can be shown to be second 
order (i.e. a simple discontinuity in specific heat) (Domb 1949); it is 
analogous to the transition in a liquid-vapour mixture at constant volume. 
For «, =} the transition takes place at C, the specific heat singularity being 
more complex. For a two-dimensional lattice it is a logarithmic infinity 
of the Onsager type. Prine 


R2 
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3.7.2. Order-disorder transitions 

When the interaction energy in a regular solution favours the formation of 
a superlattice, the eqns. (172), (173) remain valid, but the singularities 
diagram for an antiferromagnetic is now relevant to the interpretation. 
When the lattice can be divided into two equivalent sub-lattices fig. 39 is 
appropriate. When «,=} the path is A’B’ as before; for a,<}a different 
type of path again results. Examination of the series (154) for the ordered 
region shows that the dominant terms as ,, and y tend to zero are of the form 
(u-1y%)*; similarly the dominant terms in the series (143) for the disordered 
region are of the form (12%)'=(jy~%)8. Hence in both cases the relation 
(173) can only be satisfied if ~y-% tends to a finite non-zero value as y— 0. 
As y>1, p/(1 +) >a, as before. 


Fig. 42 
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Free Energy VS concentration for a Solid Solution (7’<7,). For concentrations 
between 0-1 and 0-9 the free energy curve is a horizontal line, the 
assembly having separated into two phases. 


Typical paths are drawn in fig. 39 for various values of a,. It is clear 
that for sufficiently small «, the path lies completely within the region 
D’F’C’B’E’ in which no long-range order exists; the path does not cross 
the line of singularities D’F’C’ and no critical effects will be present. On the 
other hand, for values of «, near to $, the path lies within the region of long- 
range order A’D’F’C’ for sufficiently low temperatures, and crosses the 
line of singularities to pass into the region D’'F’C’B’E’ as the temperature is 
raised; there is thus a singularity corresponding to the disappearance of 
long-range order. We deduce that a critical concentration, ae, exists below 
which there is no long-range order even at the absolute zero, and above 
which there is a normal Curie point. 
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For lattices like the triangular and F.C.C. the singularities diagram is 
fig. 40. Again we should find a critical concentration, ac, below which 
long-range order is absent. For «,>ac however, the Curie temperature 
increases to a maximum for a concentration, «,, <4, and then, decreases to 
zero for o;=4$. The discussion of this case is somewhat more tentative 
because of the lack of an adequate series for the ordered region. 


APPENDIX I 


§1. THEOREM OF PERRON AND FROBENIUS FOR SYMMETRIC MaTRICES 


Let V be a symmetric matrix of order ¢ with positive elements. The 
eigenvalues of V are all real, and with V we can associate a quadratic 
form which can be reduced to canonical form by a real orthogonal trans- 
formation. This quadratic form can therefore be written, as 

Ny (044% +... + oy, )® + Ag(Bya +... Ba )P+...+A,(O01 4+. . .6,2,)?, 
and the «;,8;...6; form an orthogonal matrix, where the \, are taken in, 
descending order of magnitude. We shall prove that if A,=A, (=1 for 
simplicity), and A3,...A,< 1, the coefficients of the quadratic form cannot 
all be positive. 

Consider the coefficient of x;x; (147). This is equal to 


ayers + BB; + Agyiy; +..-r,0:9; 
= (As— Vyeyj + Ag 1)8;8; +... +A, 1)8,0;, (174) 
in view of the orthogonality condition. All the coefficients in the right- 
hand side of (174) are negative. Thus if we define 


vit = +/(1—-Az) i; tees 6;* = 4/(1—A,)6,, 


these quantities may be regarded as components of ¢ vectors in (¢— 2) 
dimensional space, and the coefficient of 7; is minus the scalar product of 
the ith and jth vectors of this set. Our task is then to prove that if any 
t vectors are taken in (t—2) dimensional space the angles between pairs 
of vectors cannot all be obtuse. 

The theorem is elementary for t=4. We shall prove the case t=5 but 
the method can, clearly be extended to any #, and is due to Dr. 8. R. Tims 
(of King’s College, London). 

We can take the first vector §, to have coordinates (a, 0, 0) (a,>0) and 
we shall introduce the remaining vectors one at a time ensuring at each 
stage that all the angles are obtuse. The second vector must have its 
first component negative, equal to —a,(a,>0); we choose the second 
coordinate axis in the plane of &,, 6, so that &, is (— dg, b., 0) (b2 > 9). Now 
the third vector &, makes an, obtuse angle with both €, and & and hence 
its first two components are negative; the third coordinate axis can now 
be chosen, so that &3 is (—a@3, —)s3, Cg)(a3, 63, Cg>0). The remaining two 
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vectors —, and €; must make obtuse angles with €, § and €; and it will 
be seen that their coordinates are 

(—a4, —b4, — Cy) (—45, — Bs, — C5) (Aa; Oa; Ca, U5, 55, 5 > 9). 
Hence &, and &, must make an acute angle with each other. 

The generalization to the case when more than 2 of the A; are equal is 
immediate. The case A, = 1, A,= — 1 can be excluded since it would involve 
Tr (V’)~0 or 2 for large 7 according as r is odd or even. It can readily 
be shown that this is impossible since all the elements of V” are positive 
for all r. 


APPENDIX II 


§1. Exact CALCULATIONS OF THE PROPERTIES OF 
TWO-DIMENSIONAL LATTICES 


1.1. General Remarks 


We have seen in this chapter how the partition function of the Ising 
Model in zero magnetic field for certain two-dimensional lattices can be 
reduced to the evaluation of a double integral. The detailed reduction of 
this double integra] and the derivation of the thermodynamic properties 
of different lattices is given in the literature, but it seemed worth assembling 
the results together for reference purposes. Also a number of unifying 
features emerge when, the problems are treated collectively. 

We shall be primarily concerned with isotropic interactions, in which 
all directions in the lattice are equivalent, and we shall deal with the 
quadratic, triangular, honeycomb, and kagomé lattices. The partition 
function for the quadratic lattice (94) is expressed in terms of the integral 

°2n (2m 
raul ¥ In (A,—B,cos0—B,cosd)d6d¢, . . . (175) 
where A,=cosh?2K, B,=sinh2K. 

It will be seen that A,>2B, for all K, the equals sign occurring at the 
Curie point, and giving rise to a singularity in the integral at 0=¢=0. 
Similarly the partition functions for the other three lattices (127), (131) 
and (138) are expressed in terms of the integral 


1 20 (2n 
zi] i In [C, —_D, cos 0— D, cos — D, cos (0+ 4)]d0dd, (176) 
0 0 


where C, and D, have the following values: 


Triangular: C,=cosh® 2K +sinh? 2K, D,=sinh2K ; 
Honeycomb: C,=1+cosh® 2K, D,=sinh? 2K ; 
Kagomé : C= t[cosh® 2K + sinh® 2K D,=4(cosh 2K sinh’ 2K 

+ 2 cosh’ 2K sinh’ 2K + cosh? 2K sinh? 2K). 


+3 cosh? 2K], 
Again 0, >3D, for all K, the equals sign occurring at the Curie point, 
and giving rise to a singularity in the integral at 6=4¢=0. 
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The integrals (175), (176) cannot be reduced to tabulated form, and 
for the purposes of transformation it is better to consider the internal 
energy, which involves the partial derivatives with respect to A, or CO, 
of the above integrals. We therefore concentrate attention on the integrals: 


Quadratic lattice : 
2a P27 dé Rae OOD = ee 
Bee oan late 
eal ieee cond= anes fas): ee 
Other lattices : 


es dé dd 
(9, = C7>/D,). (178 
F(m) = 47? ale ie 71 — Cos 6 — cos d — cos ( +4)" ue Ay 


Using the standard result 
=f oy es Fh ME LEN 7g) 
27 Jo S—cosp 
we can, reduce both integrals (177), (178) to complete elliptic integrals of 
the first kind. 
For (177), we apply (179) to one of the variables, and use the formula 
(Bateman Manuscript Project 1953, pp. 308, 309) 


fe Agree 
eV (a—2)(B—2)(y —%)(x—8)] V [(«—y)(B—-$)] ey, 
(k= [1 kesintyy 28a, 
We find that 
G)= #0) Ee ean ee ES 


For (178), we combine the terms cos@ and cos(#+¢) and eer with 
respect to 0; we then apply (180) to the result and obtain 
2 1 2 
tlm) = a [4(v,4— 6r42 — 3) + 2v, }? #1; + 3(v1> — 6v, — aM 
ee Oe 22779) aetirees Sy ves ace (182) 

Hence we see that the energy and specific heat can be expressed in 
terms of complete elliptic integrals, but the free energy and entropy 
involve the integrals of these functions. 


1.2. Curie Points 


The Curie points of the lattices are readily determined by the trans- 
formations of §§ 3.4.2, 3.4.3, 3.4.4 and are as follows: 


Quadratic : Ge (a/2N) = 04142136,  kT'o/J =2-269185; 
Triangular : Ze = 1/1/38 =0-5773503, = kT',/J =3-640957 ; 
Honeycomb: z=(2— 1/3) =0-2679492, kT./J =1-518652; 
Kagomé: ae = [(2/4/3) — 1]J#2=0-3933198, — bT'o/ J = 2-143319. 


wee (183) 
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1.3. Critical Constants 


The detailed formulae for the internal energy per spin are obtained 
from (94), (127), (131), (138) by differentiation in the usual manner, and 
are as follows: 


Quadratic:  E£/J=—B,/2B,—-€1,(&)/2;  } 
Triangular: E/J = —D,/2D,—71o(n,)/2; 
Honeycomb: E/J = —D,/4D,—7,1o(n,)/4:; 
Kagomé : E/J = —D,/6D, —7,2o(m)/6. 


(184) 


Here the dot refers to differentiation with respect to K. At the Curie 
point we have pointed out that €, and y, have maxima, and their derivatives 
are consequently zero. Thus the critical values of the energy are deter- 
mined by the first terms on the right-hand side of (184), and are expressible 
in terms of surds. They are as follows: 


Quadratic: E,/J = — Bre —1-414214; 
Triangular : E,/J = —2; 

Honeycomb: £#,/J = —2/4/3= —1-154701; 
Kagomé : E,/J = — (1+ 21/3)/3 = — 1-488034. 


(185) 


To determine the critical values of the entropy we shall first need the 
critical values of Ind. For the quadratic lattice this involves the numerical 
value of the definite integral 


G{= a zl. fe In [1—4(cos@+cos¢)]d@d¢d, . . (186) 
which can be related to Catalan’s constant G 
G,= ~G¢—2m2=—0-2200508 . . . . . (187) 
7 


(where G = 1-?— 3-24 5-27-24... =0-915 965 594). 
For the other three lattices the corresponding definite integral is 


Qn (20 
G,= za] | In {1 —4[cos 0+ cos ¢ + cos (6+ ¢)]dédd, (188) 
0 0 


which was evaluated numerically by Houtappel as ~ —0-17642. Hence 
we have for the critical values of Ind: 


Quadratic: = InAg=In (ze! + ze) + 4G, = 0-929695 ; 

Triangular: In dAp=41n (2-3 + 32) + 4G, = 0-87959 ; 

Honeycomb: IndAc= dn [2¢-8 +. 2,3 + Bzg-2 + 82g + 8] + 4G, = 1-02506 ; 
Kagomé: In Ag = § In [ze + 182¢-? + 24 + 21262] + 4G, = 0-97479. 


(189) 
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The corresponding values of the entropy are given by 


So/k=InAg+ E.[kT, 
and are respectively 


0-30647, 0-33028, 0-26471, 0-28052. awe (EGU) 


1.4. Calculation of Thermodynamic Functions 
The internal energy, H, is given by (184), and substituting for the 
quadratic lattice, and using (181) we obtain 
E/J = —coth2K E _ i ke # (hy) | 
7 
ky =2/é; =2sinh 2K /cosh? 2K, ie es 
ky" =2tanh? 2K —1= + (1—k,2)'?. 
The specific heat per spin is given by 


d 
Cy /k= Keg EIA) 
ans (K coth K)?{2 4°(k,) —28(ky) —(1— ky") [42+ ky” H(k,)]} (192) 
where 


Am/2 
S(ky)= | (1 ky?sin? pay, 
0 
and we have used the relation (Edwards 1921, p. 399) 
(2 
i DO temad eee (bjt (2+b221). 02. 098) 
0 


The partition function can be obtained by integrating (177) with respect 
to A,, using (181), and we readily establish the formula 


In (A/2 cosh 2K) = = | "nfhl1 +(1—kesin2y)?]} dy... (194) 
7 J0 
Although this integral cannot be expressed in closed form Onsager (1944) 
derived rapidly convergent series expansions. With the notation 
ky’ =|ky"|, Inqy=aryi= — 7 4 (ky’) | 0 (hy), 
one or other of the following expressions is suitable: 


_mA=$ln2sinh2K)— tng, + (—12r—1in(1—-g?r), 


InA=4ln (2sinh 2K) + (2/n)@+ = >(-UF 
ye Ut (20 + L)(mi/7y) = exp [= (2r + Vata ee ribs) 
(2r + 1)? sinh? [(r + 4)a2/7,] 
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We have noted that k, is always <1; temperatures above and below 
the Curie point correspond to the same value of k, since k, is invariant 
under the inversion transformation (§3.4.3). The interna] energy for 
the other lattices can be obtained from (184), the appropriate », being 
invariant under the corresponding inversion transformation. For the 
triangular lattice v, of (182) takes on the simple value 1/vz, and for the 


honeycomb lattice 
1 1+2z 
L/v*z* = mak sane 
L—9}]\l—z 


in the notation of §3.4.3. The specific heat can be obtained by differen- 
tiation, and the partition function by integration of (178) with respect 
to C,. 

It is interesting to note that the elliptic integrals in (181) and (182) 
can, be reduced to forms in which the modulus is identical for quadratic, 
triangular and honeycomb lattices. Writing K,,,, for the inverse tempera- 


inv 


ture given by (43), (58) and (61) respectively, this modulus is 


sinh 2K; pity 
= = —— . [v=exp(—2K;,,)]. + 
Kt aK sate eee eee (196) 


It was this quantity which entered similarly for the spontaneous magneti- 
zation, of the three lattices. The reduction is achieved by the use of the 
Landen transformation (Bateman Manuscript project, 1953, p. 319) 


H [2](e Wl? + U2) ] = (Ly) H(i) (<1)... (197) 


For the quadratic lattice relation (43) enables us to express functions 
of K, in terms of x,. We have 


1 1 = (= 2 Kany 


1/2 
ocean “isa 
sinh2K (sinh 2K)"2 (sinh 2K)"2 NG) 4 ) Ky, (198) 


Applying this to (177) we find that 


eats be ae eee ee 
and hence using (191) and (197) 


9 
Bl T= —eoth2K| 14 = (1-14) (a) J. ai ig Ba it200} 
7 


When Ky 2 1, it is clear from (197) that «,-! must replace «, in the right- 
hand side; hence in (200) %(«,) should be replaced Dying t oh (74) 
Analogous results may be obtained for the triangular and honeycomb 


lattices although the details are more involved. For the triangular lattice 
we use the relation (58) and define 


Jy = (O78 = Die A) ee oe (200) 
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In the same manner as (198) all quantities can, be expressed in terms of 
gy, and it can be shown that the modulus of the elliptic integral in (182) 
is 2/(«,/?+ «,—?). Using (197) we find that 


2 
ie seg 91 20< 
2(7)1) mA (E-229,) Ak) («ey <1) - + . (202) 
with %(«,) being replaced by x, %(«,-) for x,>1. For the honeycomb 
lattice the relation (202) is valid, but g, is defined by analogy with (61): 


91 = (cosh 2K ,.., — 1)/(cosh 2K — 1). eet 20a) 
For both of these lattices «, is defined as for the quadratic lattice, 
k,=sinh 2K, _/sinh 2K. 


inv 


The results (200), (202) can be combined into the following formulae 
which are similar to those given by Houtappel (1950) for the three lattices : 


D) 
E/J = —coth2K + = ¢,(K) 4 (x), 
7 


Cyl Ke] —2/(sinh2K)2+ = Ge - + i) (204) 
T 


For C, we have again used the relation (193). The values of «(K) and « 
are given in, table 2, and various other properties of the lattices are sum- 
marized. The variable g, has been retained so as to demonstrate the 
symmetry of the inversion transformation. For 7'>7 >, z and v can be 
interchanged, and g, becomes g,1; thus values of the energy and specific 
heat in the region above the Curie point can, be obtained with little addi- 
tional calculation from the values below the Curie point. Near the Curie 
point the specific heat is given approximately by (128) and (132). 

The spontaneous magnetization per spin of the quadratic triangular 
and honeycomb lattices is given by the formula 


etl re eee 1 ciacs ends). 7. - (208) 

For the kagomé lattice it is (Naya 1954) 
(1432221 —22)¥2 7 12829(1 +22)8(1 + B24) TH a 
ie (1 +2?) '~ “qa 4 ap me 


1.5. Antiferromagnetic Calculations 


To derive the solution for antiferromagnetism we replace K by —K 
in all the formulae. For the quadratic and honeycomb lattices these 
formulae remain unchanged. For the triangular and kagomé lattices 
however C, and D, change in a fundamental manner. There is now no 
value of K at which the integral in (176) becomes zero; there is therefore 
no singularity and the thermodynamic functions behave smoothly and 
analytically. However as K— © (i.e. 70) C, and D, are of the same order 


(Z + 1-46), 46 ‘ (+16), “6 


CC) a a (2% 7S ee 
(T= 22). +22) al 8(c% + Dai 1(e28 — D) pot te oe 
all —e28)% z/e¢(e%— Dairle@ + D) ; 
2 (6 + 4B) z) e476, (26S + 1)(-2b— 1z p89 OE <u 
(2+%—1-2)/G=46| = (@+%—1-A(2+Z—1-2) | quodAoU0H 
s (#6 —1)¢ (z= 1-2 /(a—1_a) =» p> 5 
= ee a genes. 
Ss af (+25) gj ¢—26(,46 — T) ee Sake (0</) 
ra ; ([—2-2)/¢= 46 b=(1—2-4)(1 2-2) iepnsuelsy, 
z: al 1-(48 + 1)(Z + 46) (246 — 1) (2—,-2/(@—,-a) =4 OT > 7 
2) (g26 + 1)g—25(2°6 — T) Ot Oe, Sah; 
(z—,-2)/z=6 $=(a—,_a)(z—,-2) oryerpen?y 
2/ (z26 + 1)(,26 — 1) (2—-2)/(@—;-a) =O» OT > 7 
UOTPVULIOFSUR1Y 
i? ” h UOISIOAUT 901938 


se0144e"] quiooAoUOP] pure ‘repnsueiay, ‘omerpen’y 
jo somsodoig ormeudpoueyy, JO UoeMoleD OF SiojouVIVq “Z% IQR, 
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of magnitude; when we take out the terms corresponding to the lowest 
energy state [using the relation kInA)=S,—£,/7'| we are left with an 
integral which represents a finite entropy. Hence the zero point entropies 
are given by: 

(Triangular, Wannier 1950) 


l Qa (2a 
Ske za] | In{3 + 2[c08 0 + cos + cos (0 + 4) ]} ddd = 0-32306, 
0 0 


(Kagomé, Kano and Naya 1953) 


1 (Qa (2a 
So/k= eal In {21 —4[cos 4+ cos¢ + cos (6 + $)]d0 dd = 0-50183. 
7 0/0 
(207) 


In regard to the thermodynamic functions, for the kagomé lattice 7, 
is always positive, and (182), (184) are still valid. However, for the 
triangular lattice 7, becomes negative, and an alternative reduction to 
standard form for a pair of conjugate complex roots is required instead of 
(180). The simplest method of deriving the solution for this case is to 
continue the ferromagnetic solution, (202) for 7'> 7c using the appropriate 
transformation of the elliptic integralt. Thus v given by (60) becomes 
negative when, J <0 and hence so does «,?; we therefore use the trans- 
formation (Bateman Manuscript Project 1953, p. 319) 


He (thy [hey = hy A (ky), «2 6 ee (208) 
which in our case gives 
2 
KH (ing) = — ere HTL —np?)¥?].  . . . (209) 


The general formula (204) remains valid, with x(K) and ¢,(K) as given 
in table 2. 

1.6. Asymmetric Quadratic Net 

When the interactions differ in different directions in the lattice the 

reduction of the integrals to standard form is much more involved. Onsager 

(1944) carried out this reduction for the quadratic lattice, and we quote 
his results at the Curie point: (sinh 2K sinh 2K’=1), 

Be= —(2/m)(J’ gd 2K’ coth 2K’ + J gd 2K coth 2K), eof 210) 


Wn 
—Fo[k? =mnXo= }n2-+ (J—2yy) Incot my, +2 | an (211) 
0 


(xy, = gd 2K’ =47—gd 2K). 
Near the Curie point, 
C,,|ke~ (4/2) Ko(K +K’ sinh 2K)? sinh 2K’ — K*sinh® 2K’ gd 2K 
— 2K K'— K” sinh? 2K gd 2K~). (212) 
K,~ 4 [In (47'/|7 — Te|) —In (K coth 2K + K' coth 2K’)]. 


Sen 

+I am indebted to Dr. M. E. Fisher for this suggestion; he and Dr. M. F. 
Sykes have pointed out that there is a misprint in Houtappel’s formula for 
e,(k) in the case J <0, and that Wannier’s transformation of the elliptic integral 
is in error, 


C. Domb on 


bo 
nNG 
bo 


Notation for §3 

Xm Magnetic susceptibility of chain of m units. 

Oy, Fy, 7, Pauli spin matrices. 

(a,) ay) mean fractions of spins pointing in directions 1 and 2 
respectively. 

ayy?) Copy”) Coto12”) Carn9%”) mean fractions of pairs of spins of 
various types on sites a, b. 

of , order parameter for sites, a, b. R long-range order (o¥ ,, = R*). 

o, variable for spin at lattice point 7 taking values + 1 according as 
spin 7 points in direction 1 or 2. 

K= JT Le mer, 

MN number of nearest neighbour links of lattice. 


Rectangular lattice. Interactions J, J’. Dimensions n horizontal m vertical 
(N=mn). Nearest neighbour links 1 horizontal m vertical. 


p(r) number of closed polygons which can be constructed with r 
nearest neighbour links of the quadratic lattice. 
p(r,s) number with r horizontal s veritcal links for rectangular lattice. 
v=variable for inversion transformation, z-v; 
Ve, 2c critical values. 
v(r) number of low temperature configurations with r broken links. 
* refers to dual lattice, e.g. K*, N*, W*. 
* refers to star-triangle transformation, e.g. z*, K*. 
— refers to decoration transformation, e.g. K valu. 
Matrices C;, s; are defined by 
Cilopis. jops sage oy ae ee OE 
$,|0y A080) 04 yo, |Og oe Oe. Oey 
Quantities related to exact solution of two-dimensional Ising model (H =0) 
A,...A,... eigenvalues of matrix V. A(z) partition function. 
Yiv+++ Yon defined by (92) in terms of which A; can be simply expressed. 
Xz, Ay’, partition functions for boundaries; co, o’ boundary tensions. 
F(N,z) defined by (140), f,(z) by (148), g,(z) by (146). 
y=z', A (u,y) used in formulae for an antiferromagnetic. 
T=(1—p)/(1+ ), (4%, 7) defined by (155), ¢,(7) by (157), 


(uA,/A), 3(2)= —% (W2A,, A). 
pl 


"1, Ng fractions of spins in sub-lattices 2, 1 which have overturned. 
we critical concentration in a superlattice below which long-range 
order disappears, 
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